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Virtual Compton Scattering (VCS) on the proton has been studied at Jefferson Lab using the 
exclusive photon electroproduction reaction ep — > ep*y. This paper gives a detailed account of the 
analysis which has led to the determination of the structure functions Pll — Ptt/c and Plt, and 
the electric and magnetic generalized polarizabilities (GPs) asiQ 2 ) and /3m(Q 2 ) at values of the 
four-momentum transfer squared Q 2 = 0.92 and 1.76 GeV 2 . These data, together with the results of 
VCS experiments at lower momenta, help building a coherent picture of the electric and magnetic 
GPs of the proton over the full measured Q 2 -range, and point to their non-trivial behavior. 

PACS numbers: 13.60.-r,13.60.Fz 



I. INTRODUCTION 

The nucleoli is a composite object, and understand- 
ing its structure is the subject of intensive efforts. Its 
electromagnetic structure is cleanly probed by real and 
virtual photons. Real Compton scattering (RCS) at low 
energy gives access to the nucleon polarizabilities, which 
describe how the charge, magnetization and spin densi- 
ties in the nucleon are deformed when the particle is sub- 
jected to an external quasi-static electromagnetic field. 

Virtual Compton scattering (VCS) j*N — > jN gives 
access to the generalized polarizabilities (GPs). Being de- 
pendent on the photon virtuality Q 2 , these observables 
parametrize the local polarizability response of the sys- 
tem, i.e. they give information on the density of polariza- 
tion inside the nucleon. Experimental information on the 
GPs is obtained through the reaction of exclusive photon 
electroproduction. Several dedicated experiments on the 
proton: 

ep — > ep7 (1) 

have been performed at various Q 2 and in the low-energy 
regime. This includes the near-threshold region, where 
the center-of-mass energy W of the (jp) system is be- 
low the one-pion threshold (W < (m p + m^o), where m p 
and m^o are the proton and pion masses) , and up to the 
A(1232) resonance region. Process (1) has been stud- 
ied experimentally at MAMI [l|-0| , the Thomas Jeffer- 
son National Accelerator Facility (JLab) [5, 6] and MIT- 
Bates 0,Q. 

The results of the near-threshold VCS data analysis 
of the JLab VCS experiment E93-050, i.e. the struc- 
ture functions Pll ~ Ptt/c and Plt, and the electric 
and magnetic GPs cxe and (3m at Qf= 0.92 and 1.76 
GeV 2 , have been published elsewhere Q. However, anal- 
ysis details and cross section data were not given. This 
is the aim of the present paper, which is organized as 
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follows. After recalling briefly the theoretical concepts 
in section |nl the experimental setup is described in sec- 
tion lllll Section lTVl rcports about data analysis, including 
event reconstruction, acceptance calculation and cross 
section determination. Section [V] presents the measured 
cross section, the physics results deduced from the vari- 
ous analyses, and a discussion. A short conclusion ends 
the paper in section IVI1 



II. THEORETICAL CONCEPTS AND TOOLS 
FOR EXPERIMENTS 

This section summarizes the theoretical concepts un- 
derlying the measurements of VCS at low energy: the 
GPs, the structure functions, and the principles of mea- 
surement. For details, we refer to review papers: 0-[Lll 
(theory) and fl2T - fl6j (experiments). 



A. Generalized polarizabilities 

Polarizabilities are fundamental characteristics of 
any composite system, from hadrons to atoms and 
molecules. They describe how the system re- 
sponds to an external quasi-static electromagnetic 
field. Real Compton Scattering (RCS) yields 
for the static polarizabilities of the proton 
a E (electric) = (12.1 ±0.3 stat =p 0.4 syst ) ■ 10" 4 /m 3 
Pm (magnetic) = (1.6 ±0.4 stat ± 0A syst ) ■ 10" 4 /m 3 . 
These values are much smaller than the particle's volume 
and indicate that the proton is a very stiff object, due to 
the strong binding between its constituents. 

The formalism of VCS on the nucleon was early ex- 
plored in [lj| and the concept of generalized polariz- 
abilities was first introduced in fl9j for nuclei. The nu- 
cleon case was established within a Low Energy Theorem 
(LET), first applied to VCS by P. Guichon et al. in [H- 
This development paved the way to new experimental in- 
vestigations: it became possible to explore the spatial dis- 
tribution of the nucleoli's polarizability response, which 
is in essence the physical meaning of the GPs (see e.g. 
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1211123). 

Photon electroproduction accesses VCS via the ampli- 
tude decomposition shown in Fig. [T] T epl = T BH + 

T VCSBorn + j-VCSNonBorn^ where B JJ stan d s for the 

Bcthc-Heitler process. Formally the GPs are obtained 
from the multipole decomposition of the Non-Born am- 
plitude T VCSNonBorn taken in the "static field" limit 
Ic.m. ~* 0; where q' cm is the momentum of the final 
real photon in the 7p center-of-mass (noted CM here- 
after). The GPs are functions of q c . m ., the momentum of 
the virtual photon in the CM, or equivalently the photon 
virtuality Q 2 (see Appendix [Al for more details). After 
the work of Drechsel et al. |23l. |24{ , six independent GPs 
remain at lowest order. Their standard choice is given in 
Table HI where they are indexed by the EM transitions in- 
volved in the Compton process. Since this paper mainly 
focuses on the electric and magnetic GPs, i.e. the two 
scalar ones (or spin-independent, or non spin-flip, S=0; 
see Table [T|, we recall their definition: 

a E (Q 2 ) = -P (il ' L1)0 (Q 2 )-(^y|), 

I3 M (Q 2 ) = -P(M1,M1)0 (Q2) . (eiy?) _ 

These GPs coincide in the limit Q 2 — ► with the usual 
static RCS polarizabilities cue and Pm introduced above. 



TABLE I: The standard choice for the nucleon GPs. In the 
notation of the first column, p(p') refers to the magnetic (1) 
or longitudinal (0) nature of the initial (final) photon, L(L') 
represents the angular momentum of the initial (final) photon, 
and 5* differentiates between the spin-flip (S = 1) and non 
spin-flip (S = 0) character of the transition at the nucleon 
side. The multipole notation in the second column uses the 
magnetic (M) and longitudinal (L) multipoles. The six listed 
GPs correspond to the lowest possible order in q' c , m ,, i.e. a 
dipole final transition [L 1 = 1). The third column gives the 
correspondence in the RCS limit (Q 2 — > or q c . m . — > 0). 
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B. Theoretical models and predictions 

There are a number of theoretical models which de- 
scribe and calculate the GPs of the nucleon. They 
include: heavy baryon chiral perturbation theory 
(HBChPT) [25M23, non-relativistic quark constituent 
dispersion relations [l^, HH 32I. 
341 . effective Lagrangian model 35|, 



models |20|, [2 
lincar-cr model 




Skyrme model [36( , the covariant framework of ref. [37[ , or 
more recent works regarding GPs redefinition (38j, man- 
ifestly Lorentz-invariant baryon ChPT[H], or light-front 
interpretation of GPs [22[. 

One of the main physical interests of GPs is that they 
can be sensitive in a specific way to the various physical 
degrees of freedom, e.g. the nucleon core and the meson 
cloud. Thus their knowledge can bring novel information 
about nucleon structure. The electric CP is usually pre- 
dicted to have a smooth fall-off with Q 2 . The magnetic 
GP has two contributions, of paramagnetic and diamag- 
netic origin; they nearly cancel, making the total magni- 
tude small. As will be shown in section fV CI the available 
data more or less confirm these trends. A synthesis of di- 
verse GP predictions for the proton is presented in (3fjj . 






a) Bethe-Heitler b)VCS Born 



c)VCS Non-Born 



FIG. 1: (Color online) Feynman graphs of photon electropro- 
duction. 



The Low Energy Theorem and the structure 
functions 



The LET established in [20j is a major tool for ana- 
lyzing VCS experiments. The LET describes the photon 
electroproduction cross section below the pion threshold 
in terms of GPs. The (unpolarizcd) ep — > epj cross sec- 
tion at small q' c m is written as: 



d 5 a 



d 5 c 



.BH+Born 



0-*O + 0{q' 2 m ) ■ (1) 



The 



notation 



d 5 cr stands for d 5 a / ' dk' elab dVL' elab d£lc.M. 

its solid angle and dfl lc .m. the solid angle of the outgo- 
ing photon (or proton) in the CM; <j> is a phase-space 
factor. The (BH+Born) cross section is known and cal- 
culable, with the proton electromagnetic form factors G E 
and G P M as inputs. The ty term represents the leading 
polarizability effect. It is given by: 



where k' elab is the scattered electron momentum 



* = Vi ■ (Pll - -Ptt) 
e 



LT 



(2) 



where Pll , Ptt and Plt are three unknown structure 
functions, containing the GPs. e is the usual virtual 
photon polarization parameter and V\ , v 2 are kinemat- 
ical coefficients depending on (9 c . m .,f,fl c . m .,ip) (see 
for their full definition). The incoming photon is chosen 
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to point in the z-direction. The Compton angles arc the 
polar angle c . m , of the outgoing photon in the CM and 
the azimuthal angle ip between the leptonic and hadronic 
planes; see Fig. [2] 

The expressions of the structure functions useful to the 
present analysis are summarized here: 



Pll 
p TT 

Plt 



[Pttb P in 

2m 



■ Ge(Q 2 ) ■ /?M(Q 2 ) + [PLTspin] 



(3) 



where a em is the fine structure constant. The terms in 
square brackets are the spin part of the structure func- 
tions (i.e. containing only spin GPs) and the other terms 
are the scalar parts. The important point is that the 
electric and magnetic GPs enter only in Pll and in the 
scalar part of Plt-, respectively. 

In an unpolarized experiment at fixed Q 2 and fixed 
e, such as ours, only two observables can be determined 
using the LET: Pll — Ptt/c and Plt, he. only two 
specific combinations of GPs. To further disentangle the 
GPs, one can in principle make an e-separation of Pll 
and Ptt (although difficult to achieve), and in order to 
extract all individual GPs one has to resort to double po- 
larization |4(|. Here we perform a LET, or LEX (for Low- 
energy EXpansion) analysis in the following way: the two 
structure functions Pll — Ptt/z and Plt are extracted 
by a linear fit of the difference d 5 a exp - d 5 a BH+Born , 
based on eqs.([T]) and ([2]), and assuming the validity of the 
truncation of the expansion to 0(q' 2 m ). Then, to further 
isolate the scalar part in these structure functions, i.e. to 
access c(e{Q 2 ) and Pm{Q 2 ), a model input is required, 
since the spin part is not known experimentally. 




FIG. 2: (Color online) (ep —¥ epj) kinematics; four- 
momentum vectors notation and Compton angles (6 a . m . ,<p) 
in the yp center-of-mass. 



D. The Dispersion Relations model 

The Dispersion Relation (DR) approach is the second 
tool for analyzing VCS experiments. It is of particular 
importance in our case, so we briefly review its properties 
in this section. The DR formalism was developed by 
B.Pasquini et al. [13 S3 for RCS and VCS. Contrary 
to the LET which is limited to the energy region below 
the pion threshold, the DR formalism is also valid in the 
energy region up to the A(1232) resonance -an advantage 
fully exploited in our experiment. 



The Compton tensor is parametrized through twelve 
invariant amplitudes Fi(i = 1, • • • , 12). The GPs are ex- 
pressed in terms of the non-Born part F^ B of these am- 
plitudes at the point t = —Q 2 ,v = (s — u)/4m p = 0, 
where s, t, u are the Mandclstam variables of the Comp- 
ton scattering. The F^ B amplitudes, except for two of 
them, fulfill unsubtractcd dispersion relations. When 
working in the energy region up to the A(1232), these 
s-channel integrals are considered to be saturated by the 
ttN intermediate states. In practice, the calculation uses 
the MAID pion photo- and electroproduction multipoles 
[4lj . which include both resonant and non-resonant pro- 
duction mechanisms. 

The amplitudes F\ and F§ have an unconstrained part, 
corresponding to asymptotic contributions and dispersive 
contributions beyond ttN . For F§ this part is dominated 
by the ^-channel 7r exchange; with this input, all four 
spin GPs are fixed in the model. For F\, a main feature 
is that in the limit (t = —Q 2 , v — 0) its non-Born part is 
proportional to the magnetic GP. The unconstrained part 
of F^ B is estimated by an energy-independent function 
noted A/3, and phcnomcnologically associated with the t- 
channel cr-mcson exchange. This leads to the expression: 



p M (Q 2 ) = f3 nN (Q 2 ) + mQ' z ) 



(4) 



where /3™ N is the dispersive contribution calculated using 
MAID multipoles. The A/3 term is parametrized by a 
dipole form: 



A/3 



Dexp 



jttNi 



\Q 2 =0 



2 \2 



(i + Q7A^) 



(5) 



An unconstrained part is considered also for a third am- 
plitude, F2- Since in the limit (t = —Q 2 , v — 0) the non- 
Born part of i<2 is proportional to the sum (cue + Pm), 
one finally ends with a decomposition similar to cq.([5|) 
for the electric GP itself: 



a E (Q 2 ) 
with Aa 



a- N {Q 2 ) 

\r* ex P - fx 



Aa(Q 2 ) 



\Q 2 =0 



(6) 



ITT 



The implication for experiments is that in the DR 
model the two scalar GPs are not fixed. They depend on 
the free parameters A a and (dipole masses), which 
can be fitted from data. It must be noted that the choice 
of a dipole form in eqs.© and ^ is arbitrary: A Q and 
A^ only play the role of intermediate quantities in or- 
der to extract VCS observables, with minimal model- 
dependence. These parameters are not imposed to be 
constant with Q 2 . Our experimental DR analysis con- 
sists in adjusting A a and A^ by a fit to the measured 
cross section, separately in our two Q 2 -ranges. Then, 
in each Q 2_r ange the model is entirely constrained; it 
provides all VCS observables, at a given value of Q 2 rep- 
resentative of the range: the scalar GPs as well as the 
structure functions, in particular Pll ~ Ptt A an d Plt- 
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III. THE EXPERIMENT 

The photon clectroproduction cross section is small 
and requires a high-performance equipment to be mea- 
sured with accuracy. To ensure the exclusivity of the 
reaction, one must detect at least two of the three par- 
ticles in the final state. The chosen technique is to per- 
form electron scattering at high luminosity on a dense 
proton target, and detect in coincidence the two outgo- 
ing charged particles in magnetic spectrometers of high 
resolution and small solid angle. These devices ensure a 
clean detection and a good identification of process (1). 
Section IIIII describes how the experiment was designed 
and realized using the CEBAF electron beam and the 
JLab Hall A equipment. 



A. Apparatus 

Since the Hall was in its commissioning phase at the 
time of the data taking for this experiment (1998), not 
all devices were fully operational and the minimal num- 
ber of detectors were used. However the experiment fully 
exploited the main capabilities of the accelerator and the 
basic Hall equipment: 100% duty cycle, high resolution 
spectrometers, high luminosities. The Hall A instrumen- 
tation is described extensively in ref. I42J and in several 
thesis works related to the experiment [43144 7| . Only a 
short overview is given here, and some specific details are 
given in the subsections. 

The continuous electron beam at 4 GeV energy (un- 
polarized) was sent to a 15 cm long liquid hydrogen tar- 
get. The two High Resolution Spectrometers, noted here 
HRS-E and HRS-H, were used to detect in coincidence 
an outgoing electron and proton, respectively. After exit- 
ing the target region the particles in each HRS encounter 
successively: the entrance collimator of 6 msr, the mag- 
netic system (QQDQ), and the detector package. The 
latter consisted of a set of four vertical drift chambers 
(VDC), followed by two scintillator planes SI and S2. It 
was complemented in the HRS-E by a Cerenkov detec- 
tor and a shower counter, and in the HRS-H by a focal 
plane polarimeter. The VDCs provided particle tracking 
in the focal plane. The scintillators were the main trig- 
ger elements. They provided the timing information in 
each spectrometer and allowed to form the coincidence 
trigger. 



B. Kinematical settings and data taking 

Data were taken in two different Q 2_ran g es ; near 0.9 
and 1.8 GeV 2 . The corresponding data sets are labelled 
I and II, respectively. At Q 2 = 0.9 GeV 2 dedicated data 
were taken in the region of the nucleon resonances @. 
This leads us to split data set I into two independent 
subsets, I-a and I-b, according to the IF-range. Figure [3] 



displays the various domains covered in W , or cquiva- 
lently q' c m . Data sets I-a and II have events essentially 
below the pion threshold, while data set I-b is more fo- 
cused on the A(1232) resonance region and above. For 
the analyses presented here, emphasis will be put on data 
sets I-a and II. For data set I-b, details can be found 
already in Q, in which a nucleon resonance study was 
performed up to W = 2 GeV. Here, the lowest-VF part 
of data set I-b is analyzed in terms of GPs. Table [TT1 
summarizes our notations. 

TABLE II: The various data sets of the experiment. Columns 
2 and 3 give the ranges in Q 2 and W covered by the exper- 
iment. The fixed value of Q 2 chosen in the analyses is 0.92 
GeV 2 (resp. 1.76 GeV 2 ) for data sets I-a and I-b (resp. II). 
Columns 4 and 5 give the VK-range used in the LEX and DR 
analyses. 



data 


Q 2 -range 


W-range 


VF-range 


W-range 


set 


(GeV 2 ) 


(GeV) 


LEX (GeV) 


DR (GeV) 


I-a 


[0.85, 1.15] 


[0.94, 1.25] 


[0.94, 1.07] 


[0.94, 1.25] 


I-b 


[0.85, 1.15] 


[1.00, 2.00] 




[1.00, 1.28] 


II 


[1.60, 2.10] 


[0.94, 1.25] 


[0.94, 1.07] 


[0.94, 1.25] 





2000 






i data set I-a 

, , , ii , , , , i , , , , i 


(<Q 2 > = 0.92 GeV 2 ) 

(a) 

, , , i . , . . i . , , , 




J\ 1 data set I-b 


(<Q 2 > = 0.92 GeV 2 ) 






, , rvu 




1 

| data set II 


(<Q 2 > = 1.76 GeV 2 ) 

(c) 

, . . I , . , , I , , , , 



100 200 300 400 500 600 700 800 

q; m (MeV/c) 

J i i i i i i4 



1.0 1.2 1.4 1.6 1.8 2.0 

W (GeV) 

FIG. 3: (Color online) The range in q' c , m ,, or W, covered by 
the various data sets for the ep — > ep-y events. The vertical 
lines show the upper limit applied in the analyses: the pion 
threshold (dotted line at W = 1.073 GeV) for the LEX anal- 
yses and W — 1.28 GeV (dashed line) for the DR analyses. 

W and 5c. m. are related by W = q' c . m . + y q' 2 m. + m 2 ,. 

For each of the data sets I-a and II, the HRS-E set- 
ting was kept fixed, while the HRS-H setting was varied 
in momentum and angle. In process (1) the final proton 
is emitted in the lab system inside a cone of a few de- 
grees around the direction of the virtual photon, thanks 
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to a strong CM-to-Lab Lorentz boost. Therefore with a 
limited number of settings (and in-plane spectrometers), 
one can cover most of the desired phase space, including 
the most out-of-plane angles. As an example, Fig. 0] illus- 
trates the configuration of the HRS-H settings for data 
set I-a. In addition, in the HRS-E the momentum setting 
is chosen in order to have the VCS events in the center 
of the acceptance, i.e. near 5p/p = 0%. As a result, 
the elastic peak from ep — > ep scattering may also be 
in the acceptance of this spectrometer (at higher Sp/p), 
especially when W is low, i.e. for data sets I-a and II. 
In this case, electrons elastically scattered from hydrogen 
are seen in the HRS-E single-arm events, although they 
are kinematically excluded from the true coincidences. 
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FIG. 4: (Color online) The seventeen HRS-H settings for the 
proton detection in data set I-a. Each setting is represented 
by a box in momentum and angle. The closed curves cor- 
respond to in-plane ep — > epy kinematics at fixed values of 
q'c.m.'- 45 MeV/c (inner curve) and 105 MeV/c (outer curve). 
The ep — >• ep elastic line is also drawn at a beam energy of 
4.045 GeV. 



TABLE 111: Summary of the kinematical settings for each 
data set. The nominal beam energy is Ebeam = 4.045 GeV 
(see section HV Dl for actual values). p e and 9 e are the central 
momentum of the HRS-E spectrometer and its angle w.r.t. 
the beamline (electron side). p p and 9 P are the same variables 
on the proton side, i.e. for the HRS-H spectrometer. 



data 


Pe 


e 


p P 


P 


set 


(GeV/c) 


(deg) 


(GeV/c) 


(deg) 


I-a 


3.43 


15.4 


[0.93, 1.19] 


[45, 53] 


I-b 


[3.03, 3.26] 


[15.7, 16.4] 


[1.31, 1.53] 


[37, 45] 


II 


2.93 


23.0 


[1.42, 1.76] 


[37, 42] 



Data acquisition was performed with the CODA soft- 
ware developed by CEBAF |4^ . The trigger setup in- 
cludes several types, among which Tl and T3 are single- 
arm HRS-E and HRS-H good triggers. The T5 triggers, 
formed by the coincidence between Tl and T3, are the 
main ones used in the physics analysis. For each event the 
raw information from the detectors and the beam posi- 
tion devices is written on file. Scalers containing trigger 
rates and integrated beam charge are inserted periodi- 
cally in the datastream, as well as various parameters 
from the EPICS slow control system. Special runs were 
recorded to study spectrometer optics. 



IV. DATA ANALYSIS 

This section describes the main steps that were nec- 
essary to reach the accurate measurement of the (ep — > 
epj) differential cross section: raw-level processing, event 
reconstruction, analysis cuts, and acceptance calculation. 



A. Beam charge, target density and luminosity 

The electron beam current is measured by two reso- 
nant cavities (Beam Charge Monitors BCM) placed up- 
stream of the Hall A target. The signal of each cavity 
is sent to different electronic chains. In the experiment, 
the main measurement of the beam charge used the up- 
stream cavity and the chain consisting in an RMS-to-DC 
converter followed by a Voltage-to-Frequency converter 
(VtoF), generating pulses that are counted by a scaler. 
The content of the VtoF scaler was written on the runfile 
every 10 seconds. At the end of each run one obtained 
in this way the integrated charge of the run. The BCM 
were calibrated twice a day against the Unser monitor, 
located between the two cavities and measuring the beam 
current in absolute. The procedure also implied the of- 
fline calibration of the VtoF converter. Beam currents 
ranged from 30 to 100 fiA with an average of 70 fiA, and 
the integrated charge per run was determined with an 
accuracy of 0.5%. 

The experiment used the 15 cm long liquid hydrogen 
cell ( "beer can" ) . It was kept at a constant temperature 
T = 19.0 K and pressure P = 1.725 bar, yielding a den- 
sity po = 0.0723 g.cm -3 at zero beam current [49] . The 
beam was rastered on the target in both transverse di- 
rections in order to avoid local boiling of the hydrogen. 
Studies based on the data of this experiment [45( showed 
that density losses reached at maximum 1% for a beam 
current of 100 fiA, so the target boiling effect was con- 
sidered to be negligible and the density was taken equal 
to po m the analyses. 

The luminosity C needed for cross section measure- 
ments is obtained on a run-per-run basis. Based on the 
above considerations, it is determined with an accuracy 
of ~ ± 1%. Typical values of instantaneous luminosities 
are of the order of 2 to 4xl0 38 cm -2 -a -1 . 
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B. Rate corrections 



Event reconstruction 



The raw event rate is obtained by counting the number 
of T5 events, i.e. the coincidence triggers between the 
two spectrometers. Several correction factors have to be 
applied to this rate. 

The first correction is due to trigger inefficiency, com- 
ing from the scintillators of the detector package. It is 
obtained by studying the single-arm "junk triggers" T2 
(electron side) and T4 (hadron side), which record all 
configurations other than normal in the scintillators. The 
normal configuration (T3 or T5) is a coincidence between 
paddles in the SI and S2 planes in an allowed geometri- 
cal configuration ( "S-ray" ) , each paddle signal requiring 
the coincidence between its left and right phototubes. 
Among the junk triggers, there are some good events, 
typically with a hit missing in the scintillator paddles. 
We identify them by a "clean-up" procedure, consisting 
in the additional requirement of a valid track in the VDCs 
and a Cerenkov signal in the electron arm. The scintil- 
lator inefficiency is then defined as the number of such 
good T2 or T4 events, relative to the number of (T1+T5) 
or (T3+T5) events in the same clean-up conditions. The 
inefficiency is calculated independently for both planes 
SI and S2 in each arm. It is binned in the x (dispersive) 
and y (transverse) coordinates in each plane, to account 
for local variations. The observed inefficiency was usu- 
ally of the order of one percent, reaching occasionally 
10% locally in the electron arm [5(|. This commissioning 
problem was fixed after our experiment. 

The second correction concerns the acquisition dead- 
time. For each run, a scaler counts the number S5 of T5 
events at the output of the trigger logic. Among these 
events, only N 5 (< S 5 ) are written on file due to the 
deadtime of the [acquisition+computer] system. The cor- 
rection consists in multiplying the event rate by the ratio 
S5/N5. The deadtime depends on the beam current; it 
varied between 5% and 40% in the experiment. 

The third correction comes from the deadtime of the 
trigger electronics itself (EDT). It was not measured di- 
rectly during the experiment but determined afterwards. 
The EDT estimation is based on a fit to the actual rates 
in the scintillator paddles, using the strobe rate of each 
spectrometer. This fit was established in later experi- 
ments when the strobe rate was inserted in the datas- 
The resulting correction is of the order of 
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1-3% in our case. 

The tracking inefficiency is considered to be negligible, 
in the sense that for a real particle, the tracking algorithm 
basically always finds at least one track in the focal plane, 
which allows to process the event further. This is due to 
the good efficiency per wire in the VDCs. 

Finally, another small correction of the order of 1% is 
applied to account for the losses of protons by nuclear 
interactions in the target and spectrometer windows. 

The uncertainty on the event rate, after having cor- 
rected for all the above inefficiencies, is estimated to be 
smaller than 0.5% in relative. 



The Hall A analyzer ESPACE processes the raw detec- 
tor signals and builds all the first-level variables of the 
event: coincidence time, beam position on target, three- 
momentum vector of the detected particles at the vertex 
point, etc. Figure [5] shows a typical coincidence time 
spectrum between the HRS-E and HRS-H. The central 
peak allows to select the true coincidences. Random coin- 
cidences under the peak are subtracted using side bands. 
In the plateau one clearly sees the 2 ns microstructure of 
the beam, corresponding to Hall A receiving one third of 
the pulses of the 1497 GHz CW beam. 




150 160 170 180 190 200 210 220 230 

TCcor (ns) 

FIG. 5: Coincidence time spectrum of data set I-a. The cen- 
tral peak is 0.5 ns wide in rms. 

In the analyses presented here, particle identification 
in the detectors is basically not needed. This is because 
the kinematical settings of VCS near threshold (cf. Ta- 
ble IIIip are close to ep elastic scattering, therefore the 
true coincidences between the two spectrometers are es- 
sentially (e,p) ones. The other true coincidences that 
could be considered are of the type (e _ ,7r + ) or (ir~,p), 
coming from single or multiple pion production processes. 
However, such events either do not match the accep- 
tance settings (case of single charged-pion electropro- 
duction) or they yield missing masses which are beyond 
one pion mass, i.e. far from the VCS region of interest 
(case of multiple pion production). Therefore detectors 
such as the gas Cerenkov counter or the electromagnetic 
calorimeter in the HRS-E were essentially not used, and 
only the information from the VDCs and the scintillators 
in both arms were treated. As a verification, however, the 
analysis of data set II was performed with and without 
requiring a signal in the Cerenkov counter in the HRS-E, 
and the results were found unchanged. 

Due to the extended size of the target and the rather 
large raster size (~ ± 3-5 mm in both directions), it is 
important to know the interaction point inside the hy- 
drogen volume for each event. This point is character- 
ized by its coordinates (x v ,y v , z v ) in the Hall A labora- 
tory frame. The coordinates transverse to the beam axis, 
horizontal x v and vertical y v , are obtained from the two 
BPMs located upstream the target. It turned out that 
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for a large fraction of the data taking, the BPM infor- 
mation was accidentally desynchronized from the event 
recorded by CODA. A special re-synchronization proce- 
dure [4(1 [52| was established offline by coupling the BPM 
to the raster information (which is always synchronized 
with the event). Then the BPMs could be used, yielding 
x v and y v in absolute to better than ± 0.5 mm event per 
event. 

The calculation of the longitudinal coordinate z v re- 
quires information from the spectrometers. It is obtained 
by intersecting the beam direction with the track of one 
of the two detected particles. For this task the HRS-E 
was chosen, since it has the best resolution in horizontal 
coordinate, i.e. the variable noted yt g in the spectrometer 
frame. The resolution in y tg is excellent, about 0.6 mm 
in rms for the HRS-E (and twice larger for the HRS-H). 

The particle reconstruction proceeds as follows. In 
each arm the particle's trajectory is given by the VDCs in 
the focal plane. This "golden track" is characterized by 
four independent variables (x,y,0,(p)f p . These variables 
are combined with the optic tensor of the spectrometer 
to yield four independent variables at the target: the rel- 
ative momentum Sp/p, the horizontal coordinate yt g , and 
the projected horizontal and vertical angles, 9t g and <pt g , 
in the spectrometer frame. A fifth variable x tg character- 
izing the vertical extension of the track is calculated using 
the beam vertical position, and allows to compute small 
extended-target corrections to the dispersive variables 9 tg 
and Sp/p. The total energy of the particle is then deter- 
mined from its momentum and its assumed nature (e in 
HRS-E or p in HRS-H). It is further corrected for energy 
losses in all materials, from the interaction point to the 
spectrometer entrance. At this level the four-momentum 
vectors of the incoming electron, scattered electron and 
outgoing proton: k, k! and p', respectively, are known at 
the vertex point. One can then compute the full kinemat- 
ics of the reaction ep — > epX and a number of second- 
level variables. 

For the physics analyses, the main reconstructed vari- 
ables are Q 2 ,W, e, q c . m , on the lepton side, and the four- 
momentum vector of the missing system X: px = 
k + p — k' — p'. This four- vector is transformed from 
the laboratory frame to the CM, where one calculates 
the angles of the missing momentum vector pxc.m. w.r.t. 
the virtual photon momentum vector q c . m .. These angles, 
polar 9 c. m . and azimuthal ip (see Fig. [5]), and the modu- 
lus of the missing momentum, which represents q' c m in 
the case of VCS, are the three variables used to bin the 
cross section. 

Other second-level variables arc important for the 
event selection as well as for the experimental calibra- 
tion. The first one is the missing mass squared M x = 
(k + p — k' — p') 2 . In the experimental M x spectrum, 
a photon peak and a ir° peak are observed, correspond- 
ing to the physical processes ep — > epj and ep —> epn 
(see Fig. [7]). A cut in M\ is thus necessary to select the 
reaction channel. 

Two other variables, of geometrical nature, have 



proven to be useful. The first one, Xdif, compares two 
independent determinations of the horizontal coordinate 
of the vertex point in the Hall A laboratory frame: the 
one measured by the BPMs (x v ) and the one obtained by 
intersecting the two tracks measured by the spectrome- 
ters, and called x-2. arms . The distribution of the difference 
Xdi f = X2arms — x v is expected to be a narrow peak cen- 
tered on zero. The second geometrical variable, Y^f, will 
be described in section HV El 



D. Experimental calibration 

A lot of experimental parameters have to be well cali- 
brated. At the time of the experiment, the existing optic 
tensors of the spectrometers were not fully adapted to 
an extended target; it was necessary to optimize them. 
Using dedicated runs, new optic tensors were determined 
for the VCS analysis [HI]. They were obtained in both 
arms for our designed momentum range, and they clearly 
improved the resolution of the event reconstruction. 

A number of offsets, of either geometrical or kincmati- 
cal origin, also had to be adjusted. Among the geometri- 
cal offsets, some were given by the CEBAF surveys, such 
as the target and collimator positions. Others, such as 
the mispointing of the spectrometers, were recorded in 
the datastrcam, but their reading was not always reli- 
able and some of them had to be adjusted by software. 
One should note that not all offsets have to be known 
in absolute; what is needed is the relative coherence be- 
tween target position, beam position and spectrometer 
mispointing. The consistency checks were made on the 
distribution of the z v and Xdi / variables (defined in sec- 
tion |IVC| for real events. The main geometrical offset 
was found to be a horizontal mispointing of the HRS-E 
by 4 mm. 

The kinematical offsets consist in small systematic 
shifts in the reconstructed momentum and angles of the 
particles at the vertex point. They are mainly due to 

i) a beam energy uncertainty (the beam energy mea- 
surements described in [42j were not yet operational), 

ii) residual biases in the optic tensors, Hi) field repro- 
ducibility in the spectrometer magnets. The adjustment 
of these offsets is based on the optimization of the peaks 
in missing mass squared, in width and position, for the 
two reactions ep —¥ epj and ep — > epir° simultaneously. 
This procedure yields a coherent set of offsets for each 
setting [54[. An overview of the results is presented in 
Table [TVl All kinematical offsets were found to be small, 
except for the beam energy which was significantly below 
the nominal value from the accelerator, by about 10-16 
MeV (see Fig.©. 

E. Analysis cuts 

The offsets described above were established using 
clean event samples. However, the raw coincidences are 
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FIG. 6: The fitted offset in beam energy, AEi, eam , versus the 
setting number (time-ordered). There is one point per setting. 
The various data sets are delimited by the vertical lines. The 
horizontal line at AEt ea m = corresponds to the nominal 
beam energy from the accelerator, Ef, ea m = 4.045 GeV. 



TABLE IV: Global results for the fitted offsets on the seven 
variables: beam energy, particle momenta (p e ,p p ) and parti- 
cle angles. (ptg(e),tpt g (p) (resp. 0tg(e),0tg(p)) are the horizontal 
(resp. vertical) angles of the particle's track in the spectrom- 
eter frame. Some offsets have to be fixed in order to ensure 
the fit stability The range in brackets indicates setting- 
to-setting variations of the offsets. 



variable 


range found 
for the offset 


estimated uncertainty 
on the offset 




[-16, -10] MeV 


± 2 MeV 


Pe 


(fixed) 


± 0.3 MeV/c 


Pp 


[-1.5, +1.5] MeV/c 


± 0.5 MeV/c 


<Ptg(e) 


[0, +0.1] mr 


± 0.3 mr 


<Ptg(p) 


[-1.7, -0.7] mr 


± 0.3 mr 


0tg(e) 


[-1.6, -0.5] mr 


± 0.5 Hu- 


0tg(p) 


(fixed) 


rt 0.5 mr 



not so clean, as can be seen e.g. from the spectrum in 
the insert of Fig. [7] (top plot). The photon peak is con- 
taminated by a large broad bump centered near -15000 
MeV 2 . These events are mostly due to ep elastic scatter- 
ing with the final proton "punching through" the HRS-H 
entrance collimator. They require specific cuts in order 
to be eliminated. A key condition for the VCS analysis 
is to obtain a well-isolated photon peak in the M\ spec- 
trum ( Fig 13 histogram 5). The cuts necessary to reach 
this goal arc described below. 

First, standard acceptance cuts are applied in each 
arm. They use essentially the Hall A R- functions [55| . 
which arc a way to handle complex cuts in a multi- 
dimensional space. In the R-function approach, the prob- 
lem is transformed into the calculation, for each detected 
particle, of its "distances" to the acceptance boundaries, 
and the combination of these distances into one single 
function. This R-function takes continuous values: posi- 
tive inside the acceptance domain, negative outside, and 
equal to zero on the boundaries; it can then be used as 
a one-dimensional cut (e.g., here we require R-function 
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FIG. 7: (Color online) A sample of data set II: the experimen- 
tal spectrum of the missing mass squared at various levels of 
cuts, added successively and labelled from 1 to 5. (a): the 
raw coincidences (1) and adding the R-function cut (2). (b): 
adding the conditions W > 0.96 GeV (3), Y dif < -0.012 m 
(4) and \xdif\ < 3 mm (5) (see the text for the description of 
the variables). The insert (c) shows histogram 1 in full scale 
abscissa. 

>0). We also use additional -and largely redundant- 
contour cuts in two dimensions among the (<5, y, 9, Lp) tg 
quadruplet, and restricted apertures in the plane of the 
entrance collimators. We note that the target endcaps, 
located at ± 75 mm on the beam axis, are not seen in 
coincidence, due to the rather large HRS-H spectrometer 
angles; so a cut in z v is not necessary. The effect of the 
standard acceptance cuts is shown in Fig. [7] (histogram 
2). Clearly they are not sufficient to fully clean the M\ 
spectrum, and supplementary cuts are necessary. 

Normally, the protons coming from ep clastic scatter- 
ing are too energetic to be in the momentum acceptance 
of the HRS-H for our chosen settings. However, some of 
these protons go through the material of the entrance col- 
limator (tungsten of 80 mm thickness) where they scatter 
and lose energy, after which they enter the acceptance 
and are detected. This problem cannot be avoided, since 
VCS near threshold is by nature close to ep elastic scat- 
tering. 
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As a result, a prominent ep elastic peak is seen in the 
W^-spectrum for true coincidences, at the raw level and 
even after having applied the standard acceptance cuts, 
cf. Fig. [51 A striking evidence for "punch-through" pro- 
tons is also provided by the inserts in Fig. [8J These plots 
show the 2D impact point of the proton in the HRS-H 
collimator plane, calculated in a particular way. Here 
we do not use the information from the HRS-H, we only 
use the HRS-E informations and the two-body ep clastic 
kinematics. Knowing the vertex point (from the HRS-E 
and the beam), the beam energy and the scattered elec- 
tron angles, one can deduce the point where the proton 
from ep elastic kinematics should have hit the collima- 
tor. This hit point is characterized by its coordinates x c 
(vertical) and y c (horizontal) in the HRS-H spectrome- 
ter frame. For events in the elastic peak of Fig. [5J the 
{ x c,yc) distribution (insert (r)) reproduces faithfully the 
structure of the collimator material, proving that these 
are indeed protons punching through the collimator. The 
R-function cut is able to remove part of these events, but 
keeps the punching through the upper and lower parts of 
the collimator, as shown by insert (a). Of course our 
purpose here is only illustrative; these events are not 
a concern, since they are removed by a simple cut in 
W around the elastic peak. The result of such a cut in 
terms of missing mass squared is displayed as histogram 
3 of Fig. El 

The main concern is that there are also "punch- 
through" protons at higher W , as evidenced by insert 
(b) of Fig. [8j where an image of the collimator material 
is still present. This region in W is the far radiative tail of 
the ep elastic process; in other words it is the kinematical 
region of interest for VCS, therefore one cannot use a cut 
in W. Nevertheless these "punch-through" events must 
be eliminated, because i) they are badly reconstructed 
and ii) the simulation cannot reproduce them (our sim- 
ulation, which is used to obtain the cross section, see 
sections IIVFI and IIV Gl considers only perfectly absorb- 
ing collimators). To this aim, a more elaborate cut has 
been designed which we now describe. 

For a "punch-through" proton, the variables 
(8, y, 9, <p)t g obtained directly from the HRS-H arc 
usually severely biased, due to the crossing of a thick 
collimator. Therefore they are of little use, except for 
one particular combination: y^fj^ n = ytg + D ■ <p tg , 
where D is the distance from the target center to the 
collimator. This quantity y^Mim,™ gi ves the horizontal 
impact coordinate of the proton track in the collimator 
plane, as measured directly by the HRS-H. It is an unbi- 
ased variable, even for a "punch-through" proton. This 
is because the collimator plane is where the distortion 
happened. The reconstruction of the proton trajectory, 
which is performed backwards, from focal plane to 
target, is correct down to the entrance collimator (and 
biased further down to the target). The idea is then to 
compare this quantity yf^™ n with the "elastic" coor- 
dinate y c calculated just above. For "punch-through" 
protons the two calculations turn out to be in close 
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FIG. 8: (Color online) A sample of data set II: the experimen- 
tal W spectrum for the coincidence events after the R-function 
cut. The inserts show the proton impact in the HRS-H colli- 
mator calculated "elastically" (see text). Events in inserts (a) 
and (b) correspond to the two hatched zones of the histogram: 
the ep elastic peak (W < 0.96 GeV (a)) and a typical VCS 
region (1.0 < W < 1.073 (b)). In insert (a) a sketch of the 
tungsten collimator is drawn. The upper insert (r) shows the 
ep elastic events before applying the R-function cut. 



agreement, hence the difference Y dif = y c - y^iiirT 
peaks at zero. We point out that this comparison 
can only be done for the horizontal coordinate, not in 
vertical, due to the fact that the vertical extension x tg 
is intrinsically not measured by the spectrometers. 

Fig. [9] shows the Ydi / distribution. Clean VCS events 
cover smoothly the Ydif < region (dashed histogram), 
while experimental events (solid histogram) exhibit an 
extra-peak centered on zero. This peak corresponds to 
"punch-through" protons and is most efficiently elimi- 
nated by requiring the condition Ydif < —12 mm. The 
rejected events (insert (b)) again clearly reveal the image 
of the tungsten collimator. The retained events (insert 
(a)) show a smooth distribution in the (x c , y c ) plane, well 
reproduced by the VCS simulation. The Ydif cut is defi- 
nitely efficient in isolating a clean photon peak, as shown 
by histogram 4 of Fig. [7) 
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FIG. 9: (Color online) A sample of data set II: the experi- 
mental Ydif spectrum (see text) for coincidence events sur- 
viving the AND of the three following cuts: R-function > 0, 
-5000 < M| < 5000 MeV 2 and W > 0.96 GeV (solid his- 
togram). The inserts show the "elastically calculated" proton 
impact in the HRS-H collimator (see text) for Ydif < —0.012m 
(clean events, plot (a)) and for Ydif > —0.012 m (punch- 
through protons, plot (b)). The dashed histogram corre- 
sponds to the VCS simulation with the same three cuts. 



Lastly, to obtain histogram 5 in Fig. [7] the geometrical 
variable Xdif (cf. section ITV CI and Fig. [TU]) is selected 
around the central peak, i.e. in the interval [-3,-1-3] mm, 
completing the removal of badly reconstructed events. It 
is worth noting that these two last cuts in Ydif and Xdif 
(which arc correlated but not equivalent), owe their effi- 
ciency to the excellent spectrometer intrinsic resolution 
in y tg , already emphasized in section HV CI 

After the above cuts, a small fraction of events (< 5%) 
still have more than one track in one arm or the other. 
The number of tracks is given by the VDC algorithm to- 
gether with the parameters of the "golden track" . One 
may either keep these multi-track events, or reject them 
and renormalize the rate accordingly, based on the fact 
that these are still good events, just less clean. This sec- 
ond method was chosen, except for data set II where the 
multi-track events are in very small proportion (< 0.5%). 
Finally, events arc selected in a window in M| around 
the photon peak, typically [-5000, +5000] MeV 2 , and in 
certain VF-range. The lower bound in W corresponds to 
I'c.m. = 30 MeV/c and the upper bound is imposed de- 
pending on the type of analysis, LEX or DR (cf. Fig|3j) . 
The (very few) random coincidences that remain are sub- 
tracted. After all cuts, the final event statistics for the 
analyses are about 35000 (data set I-a), 13000 (data set 



I-b) and 25000 (data set II). 



F. Monte-Carlo simulation 

The experimental acceptance is calculated by a ded- 
icated Monte-Carlo simulation which includes the ac- 
tual beam configuration, target geometry, spectrome- 
tor acceptance and resolution effects. It is described 
in detail in [5fj| and we just recall here the main fea- 
tures. The ep —¥ epr/ kinematics are generated by sam- 
pling in the five variables of the differential cross section 
d 5 a I dk' e i ab d£l' e i ab (Kly C . m . . The scattered electron momen- 
tum and angles in the laboratory frame define the virtual 
photon, then with the angles of the Compton process in 
the CM one can build the complete 3-body kinematics. 
Events are generated according to a realistic cross sec- 
tion, the (BH+Born) one, over a very large phase space. 
The emitted electron and proton are followed in the tar- 
get materials, and the event is kept if both particles go 
through the [collimator+spcctromctcr] acceptance. One 
forms the four target variables (S, y, 9, ip)t g of each track 
and implements measurement errors on these variables 
in order to reproduce the resolution observed experimen- 
tally Finally, one proceeds to the event reconstruction 
and analysis cuts in a way similar to the experiment. 

As an example, Fig. [TU] shows the distribution of the 
variables M\ and Xdif, for two VCS data sets after all 
cuts. These variables are not sensitive to the details 
of the physics; with an infinitely good resolution they 
should be delta-functions, therefore they characterize the 
resolution achieved in the experiment. The agreement 
between the experimental and simulated data is very 
good not only in the main peak but also in the tails of 
the distributions, which is of importance as far as cuts are 
concerned. The excellent resolution achieved in missing 
mass squared allows to cleanly separate the (ep —> epj) 
and (ep — > epir ) channels. The residual contamination 
of 7T° events under the 7 peak is negligible for the settings 
analyzed here: simulation studies show that it is smaller 
than 0.5 %. 

The radiative corrections are performed along the lines 
of rcf. [57| . based on the exponentiation method. The 
simulation takes into account the internal and external 
bremsstrahlung of the electrons, because the associated 
correction depends on the acceptance and the analysis 
cuts. This allows the simulation to produce a realistic 
radiative tail in the M\ spectrum, visible on the right 
side of the peak in Fig. [TTJ-lcft. The remaining part of 
the radiative effects, due to virtual corrections plus real 
corrections which do not depend on experimental cuts, is 
calculated analytically. It is found to be almost constant 
for the kinematics of the experiment (58j : F ra d — 0.93, 
therefore it is applied as a single numerical factor, such 
that da correcte d = da raw x F ra d in each physics bin. The 
estimated uncertainty on F ra d is of the order of ±0.02, 
i.e. it induces a ±2% uncertainty on the cross section, 
globally on each point and with the same sign. 
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FIG. 10: (Color online) Data sets I-a (top) and II (bottom) 
after all cuts: comparison of experiment (solid histogram) 
and simulation (dotted histogram). (a),(c): the missing mass 
squared in the VCS region; the peak FWHM is about 1650 
MeV 2 . (b),(d): the geometrical variable Xdif (see text); the 
peak FWHM is about 1.9 mm. 



direction of q c , m . does not play a privileged role. Due 
to the narrow proton cone in the laboratory, the angular 
acceptance in the CM is almost complete for data sets 
I-a and II. This is illustrated in Fig. [TT] We note that 
the two peaks of the BH cross section, located in-plane, 
are out of the acceptance (see also Fig. [T2")) . This is on 
purpose, since in these peaks the polarizability effect in 
the cross section vanishes. When W increases, the ac- 
ceptance reduces to more backward scattering angles . 

Table [X] in Appendix [A] summarizes the bin sizes and 
the chosen fixed points in phase space. As a consequence, 
the results of the experiment are obtained at two fixed 
values of Q 2 : 0.92 and 1.76 GeV 2 . 




G. Cross section determination 

We first explain the principle of the cross section deter- 
mination in a bin, and then the chosen binning in phase 
space. In a given bin, after all cuts and corrections to 
the event rate, the analysis yields a number of experi- 
mental VCS events N exp corresponding to a luminosity 
£ e xp- Similarly, the simulation described in section TlV Fl 
yields a number of events N S i m corresponding to a lu- 
minosity Csim- The experimental cross section is then 
obtained by: 

jbEXP N exp C s j m ,5 , „s /-a 

a a- = — • — — • a cr slm {P ) , (7) 

where the factor [L s i m -d <J s i m (Po)/N s i m ] can be seen 
as an effective solid angle, or acceptance, computed by 
the Monte-Carlo method. d 5 a S i m (Po) is the cross sec- 
tion used in the simulation, at a fixed point Pq that can 
be chosen freely. As explained in [56|], this method is 
justified when the shape of the cross section d 5 a S i m is 
realistic enough, and it gives rise to a measured cross sec- 
tion (d 5 a /dk' elab dQ' elab d£l~f C . m .) at some well-defined fixed 
points in phase space. 

These points are defined by five independent variables. 
The most convenient choice w.r.t. the LET formulation 
is the set (q c . m .i q' c . m ., e , @c.m.i f)- We will work at fixed 
Qc.m. and fixed e, and make bins in the other three vari- 
ables. For the subsequent analyses, instead of the stan- 
dard (6 c .m., <P) angles, another convention (0' c m , <p' c . m .) is 
chosen. It is deduced from the standard one by a simple 
rotation: the polar axis for 9' c m is chosen perpendicu- 
lar to the lepton plane, instead of being aligned with the 
Qc.m. vector (sec Appendix [X] for more details). This new 
system of axis allows an angular binning in which the 
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FIG. 11: Accepted phase space in (9' c .m.,(p'c.m.) f° r data set 
I-a. The two crosses denote the position of the BH peaks 
and the horizontal line corresponds to in-plane kinematics 
(0' c . m . = 90° or ip = 0° and 180°). 

In eq.© the cross section d 5 a EXP is first calculated 
using the BH+Born cross section for d 5 a S i m , i.e. no po- 
larizability effect is included in the simulation. Then, 
to improve the accuracy, we include a Non-Born term 
in d 5 cr s i m , based on what we find for the polarizabilitics 
at the previous iteration. Below the pion threshold this 
Non-Born term is the first-order LEX term of eq. |T]) . For 
the region above the pion threshold, this Non-Born term 
is computed using the dispersion relation formalism, and 
the iterations are made on the free parameters of the 
model. In all cases this iterative procedure shows good 
convergence. 

H. Sources of systematic errors 

The systematic errors on the cross section come from 
three main sources: 1) overall absolute normalization, 2) 
beam energy, 3) horizontal angles of the detected parti- 
cles. 

The uncertainty in the absolute normalization has 
principally three origins: the radiative corrections known 
to ±2%, the experimental luminosity known to ±1%, 
and the detector efficiency corrections known to ±0.5% 
(see previous sections). Added in quadrature, they give 
a overall normalization error of ±2.3%, applying to all 
cross-section points with the same sign. 
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The uncertainty in the beam energy, deduced from the 
offsets study of section IIVD1 is taken equal to ±2 MeV. 
The uncertainty in horizontal angles essentially reflects 
the accuracy of the optic tensors and is taken equal to 
±0.5 mr in each arm. To study the systematic error in- 
duced by the beam energy or the horizontal angles, the 
experimental events are re-analyzed with these param- 
eters changed, one by one separately, by one standard 
deviation. One obtains in each case a set of modified 
cross-section data; in certain cases we observe a change 
of shape of the cross section. One can summarize these ef- 
fects by saying that error sources 2) and 3) taken together 
are equivalent to an average systematic uncertainty of ± 
6% (resp. ± 7%) on the cross section, for data set I-a 
(resp. II). These errors include substantial point-to-point 
correlations. 

Systematic errors on the physics observables will be 
discussed in sections IV B II and IV B 21 



I. Choice of proton form factors 

The proton clastic form factors G E and G V M are an 
important input in an analysis of VCS at low energy. 
Indeed they are needed to calculate the BH+Born cross 
section, which is at the basis of the low-energy expansion. 
Throughout these analyses the form factor parametriza- 
tion of Brash et al. [59| was chosen. It provided the first 
fit consistent with the observed de par ture from one of the 
ratio fi p G p E /G p AI in our-Q 2 range foil loll]. 

The VCS structure functions and GPs are always ex- 
tracted by measuring a deviation from the BH+Born pro- 
cess -either analytically as in the LEX approach, or in a 
more complex way as in the DR approach. This state- 
ment means that the GP extraction is sensitive to both 
cross sections, d 5 a EXP and af a BH + Born - a \% change 
on d 5 a EXP has the same impact as a 1% change on 
d 5 a BH+Born . This last cross section is not known with 
an infinite accuracy, due to uncertainties on the proton 
form factors. Therefore a systematic error should be at- 
tached to our calculation of d 5 a BH + Born , To treat it in 
a simplified way, we consider that form factor uncertain- 
ties are equivalent to a global scale uncertainty of ± 2% 
on d 5 <j BH + Born . Then, when dealing with the extrac- 
tion of the physics observables (sections IV Bl and IV CP , 
this effect can be put instead on d 5 a EXP , i.e. it can be 
absorbed in the overall normalization uncertainty of the 
experiment. Consequently, in sections IV Bl and IV CI we 
will simply enlarge the systematic error due to normal- 
ization (source # 1 in section HV Hp from ± 2.3% to ± 
3% (= quadratic sum of 2.3% and 2%). 

V. RESULTS AND DISCUSSION 

We first present the results for the photon electropro- 
duction cross section. Then the VCS structure functions 
and the GPs are presented and discussed. The main re- 



sults for these observables are contained in Tables fVl I VIII 
and US 



A. The ep — > ep7 cross section 

The experiment described here provides a unique set 
of data for VCS studies, combining altogether a large an- 
gular phase space (including out-of-planc angles) , a large 
domain in CM energy (from the threshold to the Delta 
resonance) and an access to the high-Q 2 region. Our 
cross-section data are reported in Tables IXIIII to IXVIII of 
Appendix [UJ 



1. Angular and energy dependence 

Selected samples of our results are presented in Figs. 
[12] to EI] Figure [T2l shows the measured cross section for 
the highest value of q' c m below the pion threshold (105 
McV/c). The in-plane cross section (0' cm = 90°) rises 
by seven orders of magnitude in the vicinity of the BH 
peaks, which are indicated by the two arrows. The out- 
of-plane cross section has a much smoother variation. As 
expected, the measured values exhibit a slight departure 
from the BH+Born calculation, due to the polarizabili- 
ties. The magnitude of this effect is best seen in Fig. [13] 
which depicts the deviation of the measured cross section 
relative to BH+Born: in-plane this ratio varies between 
-5% and +20%, except in the dip near ip' c _ m _ = —200° 
(or +160°) where it reaches larger values. This complex 
pattern is due to the VCS-BH interference. Out-of-planc 
the polarizability effect is much more uniform, with an 
average value of ~ —10%. 

Another selected sample of results is displayed in 
Fig.QjJ this time above the pion threshold (at q' c m =215 
MeV/c) and for backward angles of the outgoing photon. 
There, the first-order term of the LET becomes clearly 
insufficient to explain the observed cross section, while 
the calculation of the DR model, which includes all or- 
ders, performs quite well. The energy dependence of the 
cross section, i.e. the dependence in q' c m or W, is gov- 
erned by a strong rise when q' c m tends to zero due to the 
vicinity of the ep elastic scattering, and a resonant struc- 
ture in the region of the Dclta(1232). These features can 
be seen in Fig. [T51 



2. Overall normalization test 

The effect of the GPs in the photon clcctroproduction 
cross section roughly scales with the outgoing photon en- 
ergy q' c m . Therefore the physics results are determined 
essentially from the bins at high q' c m , which have the 
highest sensitivity to the GPs. At our lowest q' c m of 
45 McV/c, this sensitivity is much reduced, and one can 
test another aspect of the experiment, namely the overall 
normalization. 
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FIG. 12: (Color online) Data set I-a below the pion threshold, 
at q'c.m. = 105 MeV/c. The (ep —¥ ep~y) cross section is shown 
in-plane (8' cm = 90°, (a)) and out-of-plane (9' crn = 40°, 
(b)). The dotted curve is the BH+Born calculation. The solid 
curve includes the first-order GP effect calculated using our 
measured structure functions. The errors on the points are 
statistical only, as well as in the six next figures. The upper 
plot (c) shows the in-plane BH+Born cross section with a 
full-scale ordinate and the more traditional abscissa running 
between ip' c , m , = —180° and ip' c . m . = +180°. 
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FIG. 13: (Color online) The ratio (da EXP - 

da BH+Bo r ny da BH+Born for thg dgX& q{ ^ prg _ 

vious figure. The solid curve shows the first-order GP effect 
calculated using our measured structure functions. 



When q' c m tends to zero, d 5 a EXP formally tends to 
the known BH+Born cross section. This is a model- 
independent statement, best illustrated by the LEX ex- 
pansion of cq.([T]). At q' cm = 45 McV/c, the first-order 
term (q' cm -(p-^o), calculated using our measured values 
for the structure functions, is very small. It is about 2% 




FIG. 14: (Color online) The (ep — >• ep7) cross section for data 
sets I-a (plot (a)) and II (plot (b)) at q' c , m , = 215 MeV/c, 
in-plane (0' c m — 90°) as a function of ip' c , m ,. The dashed 
curve is the DR model calculation, with parameter values as 
fitted in the experiment. The dotted (resp. solid) curve is the 
BH+Born cross section (resp. plus a first-order GP effect). 
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FIG. 15: (Color online) Data set I-a. The (ep — ► ejry) cross 
section at the lowest q' c , m , of 45 MeV/c, for in-plane (a) and 
out-of-plane (b) kinematics. The solid curve is the (BH+Born 
+ first-order GP) cross section. The right plot (c) shows the 
reduced x 2 of the normalization test. 



of the BH+Born cross section, and it remains essentially 
unchanged when the structure functions are varied by one 
standard deviation. Therefore, at the lowest q' c m the 
comparison of the measured cross section da EXP with 
the cross section d 5 a calc calculated from cq.([T]): 



d 5 a calc = d 5 a BH+Born 



. IT, 



is essentially a test of the absolute normalization of the 
experiment. In practice, one allows da EXP to be renor- 
malized by a free factor F norm , and a \ 2 is minimized be- 
tween da EXP and da calc as a function of F norm . The test 
is performed on data sets I-a and II at the lowest q' c m ; 
the Xmin 1S always found for F norm in the range [0.99, 
1.01]. An example is given in Fig. [T5] To conclude, our 
cross-section data need very little renormalization, less 
than 1%. This means in particular that there is a good 
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consistency between the chosen parametrization of the 
proton form factors, and the way the radiative correc- 
tions arc applied to the experiment. 



B. The VCS structure functions 

As mentioned in section m the VCS structure func- 
tions and the CPs do not enter the ep — > e/ry cross 
section in the most straightforward way. A theoretical 
tool is needed to extract them from the experiment. The 
structure functions have been extracted by two different 
methods: the LEX analysis and the DR analysis. This 
section presents the methods, the results and a discus- 
sion. 



1. LEX analysis 

This analysis is based on the method described in sec- 
tion QTCj It is performed on the data sets I-a and II sepa- 
rately. An upper cut in W is imposed (W < (m p + to 7I .o), 
or q' c m < 126 MeV/c) to stay below the pion threshold. 
The LEX analysis uses the cross-section data of Tables 
Km] and EV] only. 

For each measured point in (q' c . m .,9' c .m.T P'c.m.)' one 
forms the quantity: 



AM = (cfV 



5_EXP 



d 5 , 



BH+Bc 



l )/(q' c . m .-4>) • (8) 



The V^o term of eq.([T]) is the extrapolation of AM to 
I'c.m. = in each bin in (0' c m ,ip' c m ). Below the pion 
threshold, our AM data do not exhibit any significant 
q' c m -dependence within error bars. An example is shown 
in Fig. [THJ The extrapolation to q' cm — is thus done 
simply by averaging AM over the points at q' c m = 45, 
75 and 105 McV/c. This is equivalent to neglecting the 
higher-order terms 0(q'^ m ) in this q' c m -range. 

Then a linear fit of the points, based on eq.((2|), 
yields the two structure functions Pll — Ptt/c and Plt- 
At fixed q c . m . and e, the coefficients v\ and «2 depend only 
on the Compton angles (9' c m , tp' c rn ). The good lever arm 
in vi and V2 is provided by the large coverage in these 
two angles. Figure [T71 represents the fit in terms of ^0/^2 
versus the ratio Vi/v2- Plt is given by the intercept and 
Pll — Ptt/z by the slope of the straight line fit. The 
rather good \ 2 ( c f- Table |V| confirms that higher-order 
terms 0(q' l ? m ) are small below the pion threshold. The 
values obtained for the two structure functions with their 
errors are reported in Table |V] 

The statistical errors are provided by the x 2 minimiza- 
tion. The fit can be performed on out-of-plane and in- 
plane data separately: the two corresponding types of 
results agree within statistical errors for data set I-a, but 
only within total errors (statistical + systematic) for data 
set II. 

For the systematic errors on the structure functions, 
one proceeds as for the cross section. The same sources 




100 200 100 200 100 200 100 200 

q ' c m (MeV/c) 

FIG. 16: (Color online) The q' c , m , -dependence of the quantity 
AM (see eq.©) in each angular bin. Example of data set I-a 
a t 9'c.m. = 40° (out-of-plane). The value of f' c , m , is written 
in each plot. The shaded band is the result of the LEX fit, 
"AM=constant" , within ±1<t error (statistical). The fit is 
performed on the three points below the pion threshold (filled 
circles). The other points (open circles) are above the pion 
threshold and do not participate to the LEX fit. The dashed 
curve shows the calculation of AM by the DR model, using 
the results of our DR fit (see section lV B 21) which is performed 
on all points (filled + open circles). 



of uncertainty are considered: 1) overall absolute nor- 
malization, 2) beam energy, 3) horizontal angles of the 
detected particles. The only difference is that the nor- 
malization error is now enlarged to ±3% to account for 
uncertainties in the proton form factors, as explained in 
section IIV II The LEX analysis is redone using several 
sets of modified cross section data, as described in sec- 
tion llV Hi The deviations of the structure functions w.r.t. 
the nominal analysis are recorded for all these cases, and 
finally added in quadrature. Detailed contributions to 
the systematic error are given in Table IXll of Appendix |B| 
A number of complementary systematic checks were per- 



16 




FIG. 17: (Color online) A graphical representation of the LEX 
fit for data sets I-a (plot (a)) and II (plot (b)). Each point in 
Vi/v2 corresponds to a different bin in (6' c , m ,, <p' c . m .)- The full 
circles correspond to out-of-plane data. The insert in plots 
(a) and (b) is a zoom on the in-plane data (triangles), all 
concentrated at small values of vi /vz ■ The straight line refers 
to the fit performed on all data points (in-plane + out-of- 
plane). 



GeV). However, the two-pion contribution is still small 
just above threshold; the upper limit in W in our analysis 
is taken at W=1.28 GeV, considering that the model 
calculation is able to describe the experimental data in 
this energy range (see, for example, Fig. 3 of ref. [fjj]). The 
cross-section data of Tables IXIIII to IXVIII are included. 
Globally, two different domains in W are involved: 1) the 
region of the A(1232) resonance for data set I-b, 2) the 
region essentially below the pion threshold, with a small 
extension above, for data sets I-a and II (cf. Fig. 
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TABLE V: The structure functions obtained by the LEX anal- 
ysis. The first error is statistical, the second one is the total 
systematic error. The reduced \ 2 °f the fit and the number 
of degrees of freedom are also given. 



This experiment, LEX Analyses 


data 


Q 2 


e Pll - Ptt/e Plt 


set 


(GeV 2 ) 


(GeV- 2 ) (GeV~ 2 ) 


I-a 


0.92 


0.95 1.77 ±0.24 ±0.70 -0.56 ± 0.12 ± 0.17 


II 


1.76 


0.88 0.54 ±0.09 ±0.20 -0.04 ± 0.05 ± 0.06 


I-a 




xLn = 1-22 for 32 d.o.f. 


II 




xLn = 1-50 for 31 d.o.f. 



formed, e.g. by changing the analysis cuts, or the phase- 
space points for the cross section, etc. The physics results 
obtained in these studies all stay within the systematic 
error bars of Table [V] 



2. DR analysis 

This analysis is based on the DR formalism introduced 
in section HID I It is applied to the three data sets I-a, I- 
b and II separately. The restriction to stay below the 
pion threshold is now removed. Strictly speaking, the 
DR formalism provides a rigorous treatment of the VCS 
amplitude only up to the two-pion threshold (W=1.21 



FIG. 18: (Color online) Data set I-b. The (ep — > epy) cross 
section at fixed cos 6 c . m . = —0.975 and six bins in tp (from 
Table IXVIIft . The curves show the DR model calculation for 
parameter values (A a , A j8 )=(0.70, 0.63) GeV (solid), (1.20, 
0.63) GeV (dashed), (0.70, 1.00) GeV (dotted). 

Some of our cross-section measurements above the pion 
threshold are depicted in Figs. HH and UHl Figure [TE\ 
clearly shows the excitation of the A(1232) resonance. 
The DR model gives a good description of the data and 
the various curves illustrate the sensitivity to the model 
parameters. The peak position is located at a lower mass 
than the A mass, a feature due to the VCS-BH interfer- 
ence. Figure [16] shows how the DR model reproduces 
the q' c m -dependence of the quantity AM introduced in 
section IV B II the flat behavior below the pion thresh- 
old is followed by a rise near the A resonance, where the 
higher-order terms become dominant. 

The DR formalism incorporates the BH±Born cross 
section, and a Non-Born part which contains the free 
parameters A Q and A^ (cf. section Hi Dp . The analysis 
method consists in fitting these two parameters by a \ 2 
minimization (called the "DR fit") which compares the 
model cross section to the measured one. The minimiza- 
tion cannot be solved analytically; the x 2 is computed 
on the nodes of a grid in (A a ,Ap) and its minimum is 
found numerically. For each data set a clear and single 
minimum is found, with a reasonable Xmin vamc - The 
fitted values for (A a , A^) corresponding to our three in- 
dependent data sets are reported in Table [VTI 
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The statistical errors on (A a ,Ap) are given by the stan- 
dard error ellipse at + The systematic errors are 
treated exactly as in the LEX method, i.e. by finding the 
solution in (A Q , A^) for modified cross-section sets, and 
summing quadratically the resulting variations w.r.t. the 
nominal analysis. Table IXIII of Appendix [B] gives these 
detailed contributions. One will note that the system- 
atic errors appear to be much smaller in the case of data 
set I-b. This may come from the different phase space 
coverage of this data set (in 6> c . m . and W), inducing a 
different sensitivity to the sources of errors. Remark- 
ably, all our fitted values of the (A a , A^) parameters lie 
in a narrow range: [0.63,0.79] GeV, indicating that the 
asymptotic part of the GPs (Aa, A/3 of eqs.((5|) and ([6|)) 
behaves roughly as a single dipole in the Q 2 -range of 1-2 
GeV 2 . 

TABLE VI: The fitted dipole mass parameters A a and Ap for 
the three independent data sets. The first and second errors 
are statistical and total systematic, respectively. The reduced 
X 2 of the fit and the number of degrees of freedom are also 
given. 



data 


A Q 


A|9 


set 


(GeV) 


(GeV) 


La 


0.741 ±0.040 ±0.175 


0.788 ±0.041 ±0.114 


I-b 


0.702 ±0.035 ±0.037 


0.632 ±0.036 ±0.023 


II 


0.774 ±0.050 ±0.149 


0.698 ±0.042 ±0.077 


La 


XL„ = 1-49 for 164 d.o.f. 


I-b 


Xlun = 1-34 for 328 d.o.f. 


II 


X 2 min = 1-31 for 151 d.o.f. 



Once we have the fitted values of A a and A^, the DR 
model is able to calculate the scalar GPs and the scalar 
part of the structure functions (defined in Eq.©) at the 
Q 2 under consideration. The full structure functions 
Pll — Ptt/z and Plt are then formed by adding the 
spin part. This last is parameter-free, since all spin GPs 
are fixed in the DR model. The complete DR calcula- 
tion is done separately for each data set using the inputs 
of Table EH The results for P LL - P TT /e and P LT are 
given in Table I VIII The results for Pll alone are given 
in Table IVIII1 where we have also reported the DR value 
of the spin part Ptt- 



TABLE VII: The VCS structure functions obtained by the 
DR analysis. The first error is statistical, the second one is 
the total systematic error. 



This experiment, DR Analysis 


data 


Q 2 


e Pll ~ Ptt A Plt 


set 


(GeV 2 ) 


(GeV" 2 ) (GeV" 2 ) 


La 


0.92 


0.95 1.70 ±0.21 ±0.89 -0.36 ± 0.10 ± 0.27 


I-b 


0.92 


0.95 1.50 ±0.18 ±0.19 -0.71 ± 0.07 ± 0.05 


II 


1.76 


0.88 0.40 ±0.05 ±0.16 -0.09 ± 0.02 ± 0.03 



TABLE VIII: Our result for the Pll structure function from 
the DR analysis. The value of Ptt is the parameter-free DR 
prediction. 



data 
set 


Q 2 

(GeV 2 ) 


Pll (exp.) 
(GeV" 2 ) 


Ptt (theory) 
(GeV" 2 ) 


La 


0.92 


1.19 ±0.21 ±0.89 


-0.485 


I-b 


0.92 


0.99 ±0.18 ±0.19 


-0.485 


II 


1.76 


0.24 ±0.05 ±0.16 


-0.142 



To obtain the statistical and systematic errors in Ta- 
bles [VTI] and [Vm the errors on (A a , Ap) are propagated 
to the structure functions, using the model calculation. 
One will note that the DR model exists in several ver- 
sions, each one using a different set of MAID multipolcs 
for pion electroproduction. Our analyses were done us- 
ing MAID 2000. With more recent multipole sets (MAID 
2003 or 2007), the (ep — > epj) cross section in the DR 
model changes by ~ 1-2% in our kinematics. Therefore 
the results presented here are not expected to change no- 
ticeably with the version update; they should stay largely 
within the quoted statistical error. 



3. Consistency between the different analyses 

The consistency between the two types of analysis, 
LEX and DR, can be tested only on Pll —Ptt /e and Plt 
(not on the GPs themselves) , because these two structure 
functions are the only direct outcome of the LEX anal- 
ysis. In Fig. [TO] we give a comprehensive view of all our 
measurements of the structure functions, for the three in- 
dependent data sets and the two analysis methods. The 
representation in the form of standard error ellipses in- 
dicates that error correlations between the two structure 
functions are larger in the LEX case than in the DR case. 

At Q 2 = 0.92 GeV 2 , the three measurements of Pll — 
Ptt I e agree very well. The agreement is less good on the 
three values of Plt, in particular between the separate 
DR extractions of Plt from data sets I-a and I-b, i.e. 
essentially below and above the pion threshold. These 
values become however compatible within total errors, 
including systematics. As a side remark, we note that 
a single DR analysis on the whole LL-rangc would be 
possible, by joining together data sets I-a and I-b, but it 
would mask these different results for Plt- At Q 2 = 1.76 
GeV 2 , the LEX and DR results are in mild agreement. 

Overall, Fig. Q1J] shows a rather good consistency be- 
tween the two types of extraction methods, LEX and DR, 
at each Q 2 . We also point out that the systematic error 
generally dominates in our physics results ( data sets I-a 
and II). This feature can be understood already at the 
cross section level, where the size of the systematic error 
(~ 7% of the cross section) is about half the size of the 
expected GP effect (10-15% of the cross section below 
the pion threshold). 
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FIG. 19: (Color online) The structure functions obtained by 
the two methods (LEX and DR) at Q 2 = 0.92 GeV 2 and 1.76 
GeV 2 . Each Q 2 corresponds to a solid-line polygon as marked. 
For each point, the inner ellipse is the contour at (Xmin + l)i 
yielding the standard deviation on each structure function in- 
dependently. The outer ellipse at (Xmin + 2-3) corresponds to 
a probability of 68% that both structure functions are inside 
the contour simultaneously. The statistical errors quoted in 
Tables [V] and IVHI are given by the boundaries of the inner 
contour. Dotted crosses give the size of the systematic error. 



At one given Q 2 , our LEX and DR results are obtained 
in most cases from non-independent, partially overlap- 
ping data sets. Therefore we do not propose any averag- 
ing of the points shown in Fig. Q1J1 at each Q 2 . Only the 
(I-b(DR)) and (I-a(DR)) results are truly independent 
and could possibly be averaged. 



4- Q -dependence of the structure functions 



Most of the theoretical models for GPs (section III B[) 
have a validity domain limited to low energies and low 
Q 2 . At Q 2 values of 1-2 GeV 2 , the only relevant con- 
frontation of experimental data is with the dispersive ap- 
proach, so we will focus on this model. ChPT will still 
be included as a reference in the lower-Q 2 region. 

Figure [2UJ shows the structure functions obtained in 
this experiment, together with the other existing mea- 
surements, and model calculations. The main strength of 
the JLab data is to have enlarged considerably the mea- 
sured Q -range, allowing to put in perspective different 
regions of four-momentum transfer, "high" and "low" . 

The experimental data follow the global trend of the 
models, i.e. a more or less continuous fall-off for Pr, L — 
Ptt/c-, and for Plt a rather flat behavior in the low- 
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FIG. 20: (Color online) The structure functions Pll — Prr/e 
(a) and P LT (b) measured at Bates @], MAMI and JLab 
[H (this experiment). The RCS point deduced from pj| is 
also included. The inserts are a zoom in the Q 2 -region of this 
experiment. Some points are slightly shifted in abscissa for 
visibility. The inner (outer) error bar on the points is statis- 
tical (total). The thin solid curve is the HBChPT calculation 
[27j |. The other curves show DR calculations [33] performed 
with various sets of parameters. Plot (a): A a — 0.7 GeV 
(thick solid), A Q = 1.79 GeV (dashed). Plot (b): =0.7 
GeV (thick solid), A/3 = 0.5 GeV (dashed), e = 0.9 is chosen 
to draw the curves for Pll — Prr/e. The dotted curve in the 
inserts is the spin part as given by the DR model: — Ptt/0.9 
(upper plot) and Plt spin (bottom plot). 



Q 2 region followed by an asymptotic trend to zero. At 
low Q 2 the data are in good agreement with HBChPT 
at 0(p 3 ) [13 (thin solid curve). The DR model does not 
give a parameter- free prediction of Pll — Ptt / e and Plt- 
To draw the DR curves in Fig. [20Jwe have fixed the dipole 
mass parameters A a and A^ of eqs.([5]) and ([6|), and fur- 
ther assumed that they are constant versus Q 2 . This is 
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a simplification, only aiming at a simple graphical repre- 
sentation; as explained in section Hi D[ the DR model has 
no such constraint intrinsically. The solid curve shows 
the DR calculation for typical parameter values obtained 
in our experiment: A Q = = 0.7 GeV. The dashed 
curve shows the DR calculation for other parameter val- 
ues, which agree better with some of the low-Q 2 data. 

As a general statement, there is no such single DR 
curve which goes well through all the data points, over 
the whole Q 2 -rangc. It means that a single dipole func- 
tion for the unconstrained parts Aa and A/3 of cqs.© 
and © is too limiting. This is especially true for the first 
structure function Pll — Ptt/c- all measurements are 
compatible with the thick solid curve (A a = 0.70 GeV), 
except near Q 2 ~ 0.3 GeV 2 (MAMI points) where an 
enhancement is observed in the data. This feature be- 
comes more pronounced when dealing with the electric 
GP, and will be further discussed in the next section. It 
should be noted that all measurements are performed at 
high e, around 0.9, except the MAMI points which are 
at e ~ 0.6. However this change in e can hardly account 
for the observed enhancement near Q 2 ~ 0.3 GeV 2 , since 
Ptt is expected to be a very small quantity (cf. Fig. [21]) . 

For the second structure function Plt, we note that at 
our highest Q 2 (1.76 GeV 2 ) the measured value is almost 
zero, within errors (especially the LEX result). It is sug- 
gestive of the limitations of the present extraction meth- 
ods, w.r.t. higher momentum transfers. Turning back to 
the low-Q 2 region, most models predict an extremum of 
Plt- This feature is more or less confirmed by experi- 
ment, but error bars are still large and can hopefully be 
reduced in the future. The global behavior of Plt es- 
sentially reflects the Q 2 -dependence of the magnetic GP 
Pmi which will be discussed in the next section. 

There are no measurements of the spin part of the 
structure functions, Ptt and Plt spin- Theoretical esti- 
mates are given in Fig. [5TJ In HBChPT the spin GPs 
have been calculated up to order 0(p 4 ) [62|, HU, i.e. one 
more order than for the scalar GPs, but the calculation 
does not show good convergence (cf. curves 2a and 2b 
in Fig. I2ip . This should be kept in mind when consid- 
ering Fig. the good agreement between the low-Q 2 
data and HBChPT at 0(p 3 ) may be accidental and not 
so well verified at next order. In Fig. [5T] the DR calcu- 
lation gives (parameter-free) spin structure functions of 
very small magnitude, decreasing rapidly with Q 2 and 
almost independent of the set of irN multipoles. This 
contribution is also drawn in Fig. [201 as the dotted lines. 
It amounts to 35-50% of the measured structure functions 
in the Q 2 -region of our experiment, but this percentage 
is much smaller at lower Q 2 . Experimental information 
on the spin GPs would be very valuable, but very lit- 
tle is available, due to the difficulty of such experiments 
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FIG. 21: (Color online) Theoretical predictions for the spin 
part of the measured structure functions. The DR curves 
labelled la and lb are calculated with the MAID 2000 and 
MAID 2003 multipoles, respectively. The HBChPT curves 2a 
and 2b are obtained at 0(p 3 ) [2?| and 0(p 4 ) [62], respectively. 
They are drawn up to an arbitrary value of Q^ al = 0.5 GeV 2 . 



C. The electric and magnetic GPs 

The data of this experiment allow the extraction of the 
electric and magnetic GPs of the proton at Q 2 = 0.92 and 
1 .76 GeV 2 , as an ultimate step of our analyses. In the DR 
formalism, these GPs are calculated in a straightforward 
way, once the (A a ,A^) parameters are known ("direct 
DR extraction"). On the other hand, in the LEX for- 
malism there is no such direct determination of the GPs. 
The spin structure functions Ptt and PLTspin have first 
to be subtracted from the measured ones, Pll — Ptt/z 
and Plt, using a model. For this task it is most natu- 
ral to choose the DR model, especially in our Q 2 -rangc. 
Once this subtraction is done, the last step is to remove 
the Q 2 -dependence due to the electric form factor G P E in 
the scalar part (cf.eq.([3])). Our results for aE(Q 2 ) and 
Pm{Q 2 ), following this procedure, are reported in Table 
11X1 together with the results of the direct DR extraction. 

Figure [22] summarizes the existing measurements of 
&e{Q 2 ) and /3m(Q 2 ) of the proton. It is clear that this 
picture is to some extent (DR)modcl-depcndcnt; for con- 
sistency, theoretical curves are drawn only for this par- 
ticular model. Similarly to Fig. [20] a single DR curve 
cannot reproduce all the experimental data, and an en- 
hancement can be seen in the region of the MAMI points. 
In a recent paper interpreting the GPs in the light-front 
formalism [23], the electric GP is described by adding 
to the DR calculation a Gaussian contribution centered 
near Q 2 =0.3 GeV 2 . With this parametrization, one is 
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TABLE IX: The electric and magnetic GPs extracted in this 
experiment at Q 2 = 0.92 GeV 2 and 1.76 GeV 2 . The first error 
is statistical. The second one is the total systematic error, 
obtained by propagating the errors on (A Q , A#) of Table [VTT1 
(for the DR analysis) or the errors on the structure functions 
of Table [V] (for the LEX analysis). 



data 


Q 2 


a E (Q 2 ) 


Pm(Q 2 ) 


set 


(GeV 2 ) 


(KT 4 fm 3 ) 


{1Q- 4 fm 3 ) 


DR analysis 


I-a 


0.92 


1.02 ±0.18 ±0.77 


0.13 ±0.15 ±0.42 


I-b 


0.92 


0.85 ±0.15 ±0.16 


0.66 ±0.11 ±0.07 


II 


1.76 


0.52 ±0.12 ±0.35 


0.10 ±0.07 ±0.12 


LEX analysis + [spin part subtraction by DR] 


I-a 


0.92 


1.09 ±0.21 ±0.60 


0.42 ±0.18 ±0.26 


II 


1.76 


0.82 ±0.20 ±0.44 


-0.06 ±0.17 ±0.20 



able to reproduce the measured Pll ~ Ptt/^ over the 
full Q 2 -range, and the induced electric polarization in the 
nucleon is shown to extend to larger transverse distances. 
However in [22[ no clear physical origin is associated to 
this additional and intriguing structure in q,e(Q 2 )- 




Q 2 (GeV 2 ) Q 2 (GeV 2 ) 

FIG. 22: (Color online) The world data on the electric GP 
(a) and the magnetic GP (b), with statistical (inner) and to- 
tal (outer) error bar. The solid curve is the DR calculation 
drawn for typical parameter values obtained in our experi- 
ment: A a = 0.70 GeV (left) and Ap = 0.70 GeV (right). The 
short-dashed curves show the two separate contributions to 
this calculation: the pion- nucleon intermediate states (curve 
labelled l 'irN") and the [asymptotic + beyond nN] contribu- 
tion (curve labelled "asy"). The long-dashed curve is the full 
DR calculation for other parameter values: A a = 1.79 GeV 
(left) and Ap = 0.51 GeV (right). 

In the DR model, each scalar GP is the sum of two 
terms (cf. section III D[) : the dispersion integrals, satu- 
rated by the nN contribution, and the [asymptotic ± be- 
yond ttN] contribution Aa or A/3. Figure |2"21 shows these 
two contributions separately (short-dashed curves). For 



the electric GP the [asymptotic + beyond irN] term is by 
far dominant at every Q 2 . For the magnetic GP the two 
contributions are large and of opposite sign. The ttN dis- 
persive integral is of paramagnetic nature, namely via the 
formation of the A(1232) resonance. The diamagnetism 
("asy" curve) arises from the A/3 term associated with 
the exchange of the a- meson (=[7T7r]o) in the i-channcl. 
The two terms strongly cancel, leading to an overall small 
polarizability, and a more or less pronounced extremum 
in the low-Q 2 region. 

It is well known that nucleon polarizabilities are closely 
linked to the mesonic cloud, an essential ingredient of 
nucleon structure since the first ChPT calculation of 
and Pm in RCS [671 ], The mesonic cloud is also expected 
to play an important role in the GPs. The measured 
value of the mean square electric polarizability radius of 
the proton (r 2 ), of about 2 fm 2 [H, clearly indicates that 
what is probed in VCS is a large-size structure. The 
non-trivial shape observed for the electric GP over the 
full Q 2 -range calls for further understanding, and more 
measurements of the scalar GPs in the Q 2 -region of [0-1] 
GeV 2 are needed to get a clearer picture. Such measure- 
ments are underway at MAMI [681 ] . 



VI. CONCLUSION 

The JLab E93-050 experiment was one of the first- 
generation VCS experiments, dedicated in part to the 
measurement of the generalized polarizabilities of the 
proton at high momentum transfer. It was a challenging 
task to exploit the Hall A equipment in its commission- 
ing phase to study the exclusive process ep — > ep7 and 
accurately measure its small cross section. Two methods 
have been used to extract the physics observables: one 
model-independent based on the LEX, and one model- 
dependent based on the DR formalism. The results of 
the two methods show good consistency at the level of 
the VCS structure functions Pll — Ptt/c and Plt- 

The data obtained in this experiment allow in a unique 
way to put in perspective the regions of high and low 
momentum transfer. The results are an essential piece 
to build a more complete picture of the electric and mag- 
netic GPs of the proton as a function of Q 2 , i.e. ulti- 
mately the nucleon's polarization response as a function 
of the distance scale. The electric GP does not seem 
to have a smooth fall-off, and the behavior of the mag- 
netic GP quantifies the detailed contributions of para- 
and diamagnetism in the proton. Experimental data are 
still scarce, and more measurements are desirable in or- 
der to improve our understanding of these fundamental 
observables. This is especially true at low Q 2 , where the 
prominent role of the mesonic cloud can be probed. New 
VCS experiments, together with new RCS experiments 
and theoretical developments in the field, should provide 
a step forward in our understanding of the electromag- 
netic structure of the nucleon. 
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Appendix A: Kinematics and Binning 

The GPs depend on q c . m ., or cquivalently on the 
four-momentum transfer squared Q 2 taken in the limit 



@. This variable is defined as Q 2 = 2m p ■ 
r qc.m. 2 — Trip). It has been denoted Q 2 through- 
out the paper for simplicity. 

Two angular systems in the CM are described in sec- 
tion IIV Gt "standard" and "rotated" . In the standard 
system, the polar angle 6 c , m . is measured w.r.t. the z 
axis aligned with the q c . m , vector. In the rotated sys- 
tem, the polar angle 0' cm is measured w.r.t. the z' axis 
orthogonal to the leptonic plane; see Fig. [23l In-plane 
kinematics correspond to 0' c m = 90°, or to (p = 0° and 
180°. The conversion formulae between the two systems 
are the following: 



cosV cm = sinc/ c . m . • smtp , 

cos tp' cm = cos6> c . m . / sin0' CTO , 

s'mtp' cm = sin0 c . m . • costp I , &va.Q' cm 



above. For Q' c m , the up-down symmetry w.r.t. the lep- 
tonic plane allows to sum the intervals as indicated (U 
symbol). This up+down sum is also applied in the tp- 
bins of data set I-b. 



TABLE X: The phase-space variables (first column), the bin 
size and the points Pq at which the cross section is determined 
(second column). These points are chosen at the middle of 
the bin, except for the last bin in q' c .m. and the first bin in 

"cm. • 

Data sets I-a and II 
e fixed at 0.950 (0.879) for data set I-a (II) 

q c . m . fixed at 1.080 (1.625) GeV/c for data set I-a (II) 
q' c . m . 6 bins: regular from 30 to 180, + [180, 250] MeV/c 

6 points: 45, 75, 105, 135, 165, 215 MeV/c 
0' c . m . 2 bins: [0°,60°] U [120°, 180°]; [60°, 90°] U [90°, 120°] 

2 points: 40° (=out-of-plane) and 90° (=in-plane) 
ip' c . m . 20 bins: regular from -180° to +180° 

20 points: -171°, -153°, -135°, ... ,+153°, +171° 



Data set I-b (cf. ref. [6]) 



klab 

Q 2 
w 



fixed at 4.032 GeV 
fixed at 1.0 GeV 2 

15 bins: regular from 0.98 GeV to 1.28 GeV 

15 points: 0.99, 1.01, 103 , ... , 1.25, 1.27 GeV 

3 bins: [-1, -0.95], [-0.95, -0.80], [-0.80, -0.50] 

3 points: -0.975, -0.875, -0.650 

6 bins: regular from to 180° 

6 points: 15°, 45°, 75°, 105°, 135°, 165° 



(Al) 



Appendix B: Systematic errors 




9 





I M 
1 JF 












/ plane 


qcm / 



This Appendix gives details on the systematic errors on 
the primary observables extracted from the cross section 
data: the structure functions in the LEX analysis (Ta- 
ble IXI|) and the (A Q , A/3) parameters in the DR analysis 
(Tabic [XTTjl . The systematic errors on all other extracted 
observables are derived from these ones by error propa- 
gation. The origin of the 3% normalization error in these 
Tables is explained in section IIV II 



Appendix C: Cross section Tables 



FIG. 23: (Color online) The angular systems to measure the 
polar and azimuthal angles of the momentum vector q c .m. in 
the (7p) center-of-mass. 

Table |X] gives the binning information for the three 
data sets. The variables which are kept fixed arc differ- 
ent for (I-a and II) and for (I-b). The three first bins 
in q' c m are below the pion threshold, the other ones are 



This Appendix gives our measured values of the pho- 
ton electroproduction cross section. Tabic IXIIII (resp. 
IXIVp is for data set I-a (resp. II) below the pion thresh- 
old. Table IXVI (resp. IXVip is for data set I-a (resp. 
II) above the pion threshold. In these four Tables, the 
errors are statistical only (in r.m.s); the systematic error 
is discussed in the text (section IIV HI) . 

The cross-section values for data set I-b can be found 
partly in ref. Q, but there some settings of data set I-a 
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TABLE XI: Detailed systematic errors on the structure func- 
tions in the LEX analysis. To obtain the total systematic 
error, the three partial contributions are first symmetrized 
(+/— ) and then added in quadrature. 



source of 


A{P LL - P TT /e) 


A(P LT ) 


systematic error 


(GeV~ 2 ) 


(GeV~ 2 ) 


Data Set I-a 


normalization (± 3 %) 


+0.505 -0.505 


+0.046 -0.046 


beam energy (± 2 MeV) 


+0.391 -0.354 


+0.132 - 0.024 


spec. angle (± 0.5 mr) 


+0.301 -0.301 


+0.148 -0.148 


Total syst. error 


± 0.696 


± 0.174 


Data Set II 


normalization (± 3 %) 


+0.142 -0.142 


+0.004 -0.004 


beam energy (± 2 MeV) 


+0.139 -0.074 


+0.017 -0.005 


spec. angle (± 0.5 mr) 


+0.096 -0.096 


+0.054 -0.054 


Total syst. error 


± 0.202 


± 0.055 



were also included. Thus, for the sake of completeness, 
the cross-section values corresponding to the DR analy- 
sis (I-b) of the present paper are reported here in Table 
IX VII I It should be noted that this cross section has been 
determined at Q 2 = 1.0 GeV 2 (instead of 0.92 GeV 2 ), 
as part of a wider study program at Q 2 = 1 GeV 2 [a]. 
To obtain the VCS observables at Q 2 = 0.92 GeV 2 from 
the "I-b" DR fit, we have explicitcly assumed that the 
(A a , A^) parameters have no significant variation locally 
in Q 2 . Then we simply use the parameter values of Ta- 
ble ED (I-b part), fitted at Q 2 = 1.0 GeV 2 , as inputs to 
the DR calculation of the VCS observables at Q 2 = 0.92 
GeV 2 . 



TABLE XII: Detailed systematic errors on the parameters A a 
and A/3 in the DR analysis. To obtain the total systematic 
error, the three partial contributions are first symmetrized 
(+/— ) and then added in quadrature. 



source of 


A (A Q ) 


A(A /3 ) 


system. error 


(GeV) 


(GeV) 


Data Set I-a 


normalization (± 3 %) 


+0.105 -0.099 


+0.068 -0.053 


beam energy (± 2 MeV) 


+0.179 -0.067 


+0.068 - 0.018 


spec. angle (± 0.5 mr) 


+0.072 -0.072 


+0.087 -0.087 


Total syst. error 


± 0.175 


± 0.114 


Data Set II 


normalization (± 3 %) 


+0.096 -0.127 


+0.057 -0.056 


beam energy (± 2 MeV) 


+0.118 -0.047 


+0.044 -0.017 


spec. angle (± 0.5 mr) 


+0.052 -0.052 


+0.041 -0.041 


Total syst. error 


± ± 0.149 


± 0.077 


Data Set I-b 


normalization (± 3 %) 


+0.031 -0.041 


+0.018 -0.014 


beam energy (± 2 MeV) 


+0.017 -0.003 


+0.011 -0.012 


spec. angle (± 0.5 mr) 


+0.005 -0.005 


+0.012 -0.012 


Total syst. error 


± ± 0.037 


± 0.023 



Ascii files of Tables IXffll to IXVlTl will be added as aux- 
iliary files to this paper. They arc also available at URLs: 
|http://userweb.jlab.org/~helene/paper_vcs-gps-2012/all-ascii-table£ 
and http : / / clrwww. in2p3 . fr / sondcm/E93050-tables- 
GPS, or upon request to the authors. 
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TABLE XIII: The measured (ep — > epj) cross section d 5 a / dk' elab dQ! £lah dQ, lc . m . , in pb/(GeV sr 2 ), for data set I-a below the 
pion threshold. The out-of-plane (resp. in-plane) data correspond to 6' c , m , = 40° (resp. 90°). The error Act is statistical only. 
The (ep — > epy) kinematics are entirely determined by the five variables (g c .m., e, q' c . m .j "c.m.i ¥>e.ro.)- 



at fixed q c . m . = 1.080 GeV/c, fixed e = 0.950 and fixed Q 2 = 0.92 GeV 2 



@cm. 




"cm. 


tp'c.m. 


d 5 a 


± AtJstat 


d 5 a ± A<Jstat 


d 5 


c ± Aer stat 


(deg) 


(deg) 


(deg) 


(deg) 


a t 9c. m. 


— 


45 MeV/c 


a l Qc.m 




: 75 MeV/c 


Qc.m 




: 105 MeV/c 


50.6 


82.5 


40 


9 


105.5 


± 


12.0 


68.0 


± 


5.0 


48.7 


± 


3.9 


55.1 


69.1 


40 


27 


102.6 


± 


17.2 


75.5 


± 


6.1 


53.9 


± 


4.5 


63.0 


59.3 


40 


45 


141.8 


± 


31.0 


114.5 


± 


10.2 


75.0 


± 


6.2 


73.0 


53.2 


40 


63 


130.5 


± 


51.3 


123.9 


± 


16.3 


117.1 


+ 


11.4 


84.2 


50.3 


40 


81 


263.7 


± 


116.3 


238.1 


± 


33.9 


187.2 


± 


19.3 


95.8 


50.3 


40 


99 


536.3 


± 


233.1 


417.2 


± 


62.5 


243.8 


± 


28.2 


107.0 


53.2 


40 


117 


1210.1 


± 


425.1 


445.1 


± 


67.2 


300.8 


+ 


32.1 


117.0 


59.3 


40 


135 


1389.4 


± 


306.0 


491.4 


± 


59.2 


291.8 


± 


25.2 


124.9 


69.1 


40 


153 


619.3 


± 


127.0 


458.1 


± 


39.1 


247.0 


+ 


15.8 


129.4 


82.5 


40 


171 


570.7 


± 


80.2 


340.9 


± 


21.7 


219.7 


± 


10.5 


129.4 


97.5 


40 


-171 


556.9 


± 


50.4 


318.1 


± 


13.9 


204.4 


± 


8.4 


124.9 


110.9 


40 


-153 


576.2 


± 


34.1 


286.3 


± 


10.3 


191.9 


± 


9.2 


117.0 


120.7 


40 


-135 


443.0 


± 


21.2 


248.7 


± 


8.9 


149.2 


± 


8.6 


107.0 


126.8 


40 


-117 


367.7 


± 


15.9 


196.3 


± 


8.0 


130.4 


± 


8.1 


95.8 


129.7 


40 


-99 


281.8 


± 


12.7 


155.0 


± 


7.1 


93.2 


+ 


6.1 


84.2 


129.7 


40 


-81 


242.8 


± 


11.7 


119.6 


± 


6.0 


88.7 


± 


5.4 


73.0 


126.8 


40 


-63 


166.5 


± 


9.6 


106.4 


± 


5.5 


57.1 


± 


3.9 


63.0 


120.7 


40 


-45 


120.2 


± 


8.3 


73.8 


± 


4.3 


50.7 


± 


3.6 


55.1 


110.9 


40 


-27 


108.3 


± 


8.6 


60.9 


± 


4.2 


48.5 


± 


3.5 


50.6 


97.5 


40 


-9 


107.7 


± 


9.5 


55.8 


± 


4.0 


36.5 


± 


3.1 


9.0 


0.0 


90 


9 


9.5 


± 


2.4 


15.6 


± 


2.3 


23.1 


± 


3.3 


27.0 


0.0 


90 


27 


7.2 


± 


3.0 


12.8 


± 


2.2 


12.6 


± 


2.3 


153.0 


0.0 


90 


153 


33.6 


± 


8.6 


11.9 


± 


3.5 


8.8 


± 


2.7 


171.0 


0.0 


90 


171 


11.3 


± 


2.5 


10.6 


± 


1.4 


9.8 


± 


1.0 


171.0 


180.0 


90 


-171 


258.4 


± 


25.8 


128.6 


± 


7.6 


89.6 


± 


5.1 


153.0 


180.0 


90 


-153 


553.1 


± 


30.7 


286.1 


± 


10.5 


178.7 


± 


10.0 


135.0 


180.0 


90 


-135 


468.1 


± 


18.8 


266.0 


± 


10.1 


167.6 


± 


11.7 


117.0 


180.0 


90 


-117 


402.3 


± 


14.1 


231.1 


± 


9.8 


149.4 


± 


10.3 


99.0 


180.0 


90 


-99 


286.8 


± 


11.2 


159.0 


± 


8.1 


99.3 


+ 


6.3 


81.0 


180.0 


90 


-81 


194.5 


± 


9.0 


95.6 


± 


5.4 


62.6 


± 


4.0 


63.0 


180.0 


90 


-63 


121.8 


± 


7.2 


62.4 


± 


4.1 


50.8 


± 


3.6 


45.0 


180.0 


90 


-45 


74.1 


± 


5.7 


46.6 


± 


3.3 


30.2 


+ 


2.8 


27.0 


180.0 


90 


-27 


39.7 


± 


4.1 


31.4 


± 


2.6 


23.2 


± 


2.8 


9.0 


180.0 


90 


-9 


17.7 


± 


2.7 


27.5 


± 


2.8 


18.6 


± 


2.7 
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TABLE XIV: Same as previous Table but for data set II below the pion threshold. 



at fixed q c . m . = 1.625 GeV/c , fixed e = 0.879 and fixed Q 2 = 1.76 GeV 2 



$c.m. 


9 


cm. 


f'c.m. 


d a ± Aa s tat 


d E 


V ± Aa a tat 


d 


r 'a 


± A(Jstat 


(deg) 


(deg) 


(deg) 


(deg) 


a f Qc.TTl 




: 45 MeV/c 


at q c , 


= 75 MeV/c 


at q' c/ 


■'11. 


= 105 MeV/c 


50.6 


82.5 


40 


9 


9.3 


± 


1.2 


5.6 


± 0.6 


2.9 


± 


0.4 


55.1 


69.1 


40 


27 


8.0 


± 


1.5 


5.0 


± 0.7 


3.7 


± 


0.6 


63.0 


59.3 


40 


45 


6.9 


± 


2.0 


7.1 


± 1.2 


5.4 


± 


0.9 


73.0 


53.2 


40 


63 


16.5 


± 


4.8 


8.8 


± 1.7 


7.6 


± 


1.6 


84.2 


50.3 


40 


81 


10.3 


± 


5.6 


15.7 


± 3.3 


9.1 


± 


1.9 


95.8 


50.3 


40 


99 


80.2 


± 


21.1 


30.5 


± 5.4 


25.5 


± 


4.1 


107.0 


53.2 


40 


117 


87.2 


± 


23.2 


40.1 


± 6.8 


26.1 


± 


4.1 


117.0 


59.3 


40 


135 


122.5 


± 


25.9 


42.5 


± 6.0 


31.1 


± 


3.8 


124.9 


69.1 


40 


153 


140.5 


± 


25.3 


51.1 


± 6.1 


27.0 


± 


2.8 


129.4 


82.5 


40 


171 


62.2 


± 


12.0 


38.9 


± 3.7 


25.1 


± 


2.3 


129.4 


97.5 


40 


-171 


74.3 


± 


9.2 


34.3 


± 2.3 


20.1 


± 


1.5 


124.9 


110.9 


40 


-153 


51.7 


± 


4.6 


33.1 


± 1.6 


18.6 


± 


0.9 


117.0 


120.7 


40 


-135 


42.4 


± 


2.9 


25.4 


± 1.2 


15.3 


± 


0.8 


107.0 


126.8 


40 


-117 


36.1 


± 


2.2 


18.9 


± 0.9 


10.6 


± 


0.7 


95.8 


129.7 


40 


-99 


28.4 


± 


1.7 


12.8 


± 0.7 


7.0 


± 


0.6 


84.2 


129.7 


40 


-81 


19.4 


± 


1.6 


10.1 


± 0.7 


6.5 


± 


0.6 


73.0 


126.8 


40 


-63 


14.5 


± 


1.5 


7.2 


± 0.6 


3.7 


± 


0.4 


63.0 


120.7 


40 


-45 


10.6 


± 


1.2 


5.5 


± 0.5 


3.9 


± 


0.5 


55.1 


110.9 


40 


-27 


8.4 


± 


1.1 


6.3 


± 0.6 


3.2 


± 


0.5 


50.6 


97.5 


40 


-9 


9.0 


± 


1.3 


3.9 


± 0.5 


3.6 


± 


0.5 


9.0 


0.0 


90 


9 


1.5 


± 


0.5 


1.3 


± 0.4 


1.1 


± 


0.3 


27.0 


0.0 


90 


27 


1.2 


± 


0.7 


0.7 


± 0.3 


1.3 


± 


0.4 


171.0 


0.0 


90 


171 


5.3 


± 


1.2 


1.6 


± 0.4 


0.9 


± 


0.2 


171.0 


180.0 


90 


-171 


39.9 


± 


7.3 


25.5 


± 2.4 


15.8 


± 


1.1 


153.0 


180.0 


90 


-153 


64.1 


± 


6.0 


44.7 


± 2.1 


26.4 


± 


1.1 


135.0 


180.0 


90 


-135 


65.6 


± 


3.7 


32.2 


± 1.3 


19.9 


± 


0.9 


117.0 


180.0 


90 


-117 


42.6 


± 


2.1 


20.8 


± 0.9 


12.7 ± 


0.8 


99.0 


180.0 


90 


-99 


24.4 


± 


1.3 


12.8 


± 0.7 


8.3 


± 


0.7 


81.0 


180.0 


90 


-81 


14.1 


± 


0.9 


8.2 


± 0.6 


4.5 


± 


0.5 


63.0 


180.0 


90 


-63 


10.7 


± 


0.8 


5.5 


± 0.5 


3.5 


± 


0.5 


45.0 


180.0 


90 


-45 


7.4 


± 


0.8 


3.8 


± 0.4 


2.1 


± 


0.4 


27.0 


180.0 


90 


-27 


3.4 


± 


0.6 


3.2 


± 0.4 


2.1 


± 


0.4 


9.0 


180.0 


90 


-9 


2.8 


± 


0.7 


1.8 


± 0.4 


0.9 


± 


0.3 
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TABLE XV: Same as previous Table but for data set I-a above the pion threshold. 



at fixed q c . m . = 1.080 GeV/c , fixed e = 0.950 and fixed Q 2 = 0.92 GeV 2 



@c.m. 


<P 


a' 

° cm. 




d 5 


c ± A<T stat 


d 5 




d (7 ± AtJstat 


(deg) 


(deg) 


(deg) 


(deg) 


a t Qc.m 




: 135 MeV/c 


a t Qc.m 




= 165 MeV/c 


at q' c . m , = 215 MeV/c 


50.6 


82.5 


40 


9 


38.8 


± 


5.2 


32.8 


± 


20.6 




55.1 


69.1 


40 


27 


44.5 


± 


5.9 


53.6 


± 


18.1 




63.0 


59.3 


40 


45 


54.2 


± 


6.1 


43.8 


± 


10.1 




73.0 


53.2 


40 


63 


91.4 


± 


11.6 


97.5 


± 


21.9 




84.2 


50.3 


40 


81 


91.1 


+ 


15.5 


168.3 


± 


42.7 




95.8 


50.3 


40 


99 


175.5 


+ 


26.2 


143.0 


± 


38.8 


159.6 ± 95.1 


107.0 


53.2 


40 


117 


229.0 


± 


28.7 


169.4 


± 


33.4 


168.8 ± 72.4 


117.0 


59.3 


40 


135 


193.8 


+ 


19.9 


153.2 


± 


19.6 


157.1 ± 24.9 


124.9 


69.1 


40 


153 


184.4 


± 


14.0 


124.1 


± 


11.3 


119.7 ± 14.6 


129.4 


82.5 


40 


171 


178.5 


± 


12.1 


125.2 


± 


12.6 


125.5 ± 18.3 


129.4 


97.5 


40 


-171 


121.4 


+ 


9.9 


127.1 


± 


14.3 


155.6 ± 39.6 


124.9 


110.9 


40 


-153 


153.8 


± 


13.4 


81.2 


± 


19.4 




117.0 


120.7 


40 


-135 


90.2 


+ 


11.3 


141.0 


± 


29.3 




107.0 


126.8 


40 


-117 


74.7 


+ 


8.6 


80.7 


± 


28.3 




95.8 


129.7 


40 


-99 


67.3 


± 


7.0 


71.9 


± 


33.1 




84.2 


129.7 


40 


-81 


69.9 


+ 


8.2 


36.0 


± 


33.9 




73.0 


126.8 


40 


-63 


32.2 


+ 


4.7 










63.0 


120.7 


40 


-45 


45.5 


± 


6.0 










55.1 


110.9 


40 


-27 


33.9 


+ 


5.3 










50.6 


97.5 


40 


-9 


41.2 


± 


5.9 










9.0 


0.0 


90 


9 


14.9 


± 


5.6 










27.0 


0.0 


90 


27 


22.7 


± 


6.2 










45.0 


0.0 


90 


45 


18.1 


± 


8.2 










135.0 


0.0 


90 


135 














176.9 ± 75.4 


153.0 


0.0 


90 


153 


7.2 


+ 


1.7 


16.5 


± 


3.0 


80.1 ± 12.6 


171.0 


0.0 


90 


171 


14.6 


± 


1.8 


31.1 


± 


4.5 


79.7 ± 17.9 


171.0 


180.0 


90 


-171 


84.5 


± 


7.8 


94.0 


± 


11.8 


216.1 ± 70.3 


153.0 


180.0 


90 


-153 


137.6 


± 


12.8 


176.5 


± 


52.7 




135.0 


180.0 


90 


-135 


112.2 


± 


16.3 










117.0 


180.0 


90 


-117 


61.8 


+ 


9.8 










99.0 


180.0 


90 


-99 


56.3 


± 


14.4 










81.0 


180.0 


90 


-81 


45.8 


± 


20.3 










27.0 


180.0 


90 


-27 


28.1 


± 


10.4 










9.0 


180.0 


90 


-9 


30.8 


± 


9.3 
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TABLE XVI: Same as previous Table but for data set II above the pion threshold. 







at fixed q c . m . 


= 1.625 GeV/c , fixed e 


= 0.879 and fixed Q 2 = 


1.76 GeV 2 




@c.vn. 


<P 


a' 

° cm. 




d 5 a ± AtJstat 


d a ± A(T B tat 




± AtJstat 


(deg) 


(deg) 


(deg) 


(deg) 


at q' c . m . = 135 MeV/c 


at q'c.m. = 165 MeV/c 


at qc.m. 


= 215 MeV/c 


50.6 


82.5 


40 


9 


2.3 ± 0.5 








55.1 


69.1 


40 


27 


3.8 ± 0.7 


2.8 ± 1.1 






63.0 


59.3 


40 


45 


4.1 ± 0.8 


3.4 ± 1.1 






73.0 


53.2 


40 


63 


5.0 ± 1.4 


5.3 ± 2.0 






84.2 


50.3 


40 


81 


7.9 ± 2.1 


4.2 ± 2.4 






95.8 


50.3 


40 


99 


10.8 ± 2.6 


19.7 ± 5.9 


9.3 ± 


5.1 


107.0 


53.2 


40 


117 


18.3 ± 2.8 


9.1 ± 2.1 


14.4 ± 


3.4 


117.0 


59.3 


40 


135 


21.3 ± 2.7 


13.9 ± 2.3 


9.9 ± 


2.1 


124.9 


69.1 


40 


153 


18.0 ± 2.0 


15.0 ± 1.8 


12.3 ± 


1.0 


129.4 


82.5 


40 


171 


17.1 ± 1.3 


13.6 ± 1.1 


8.8 ± 


0.8 


129.4 


97.5 


40 


-171 


18.2 ± 1.0 


13.2 ± 0.9 


8.6 ± 


0.7 


124.9 


110.9 


40 


-153 


14.3 ± 0.8 


9.7 ± 0.8 


8.0 ± 


1.2 


117.0 


120.7 


40 


-135 


9.7 ± 0.7 


9.5 ± 1.1 


8.7 ± 


2.6 


107.0 


126.8 


40 


-117 


7.5 ± 0.7 


5.0 ± 1.0 






95.8 


129.7 


40 


-99 


6.2 ± 0.7 


8.0 ± 2.2 






84.2 


129.7 


40 


-81 


5.6 ± 0.8 


2.6 ± 1.2 






73.0 


126.8 


40 


-63 


4.0 ± 0.7 


4.9 ± 2.9 






63.0 


120.7 


40 


-45 


3.3 ± 0.7 


5.2 ± 3.0 






55.1 


110.9 


40 


-27 


2.0 ± 0.5 


2.3 ± 1.5 






50.6 


97.5 


40 


-9 


3.9 ± 0.9 








9.0 


0.0 


90 


9 


1.9 ± 1.0 








27.0 


0.0 


90 


27 


0.7 ± 0.4 








153.0 


0.0 


90 


153 






3.7 ± 


1.3 


171.0 


0.0 


90 


171 


1.4 ± 0.2 


1.7 ± 0.2 


4.2 ± 


0.4 


171.0 


180.0 


90 


-171 


12.7 ± 0.8 


10.1 ± 0.7 


9.5 ± 


0.7 


153.0 


180.0 


90 


-153 


19.1 ± 0.9 


14.4 ± 1.0 


13.6 ± 


1.9 


135.0 


180.0 


90 


-135 


13.1 ± 0.8 


10.3 ± 1.5 






117.0 


180.0 


90 


-117 


8.6 ± 0.9 








99.0 


180.0 


90 


-99 


6.2 ± 1.1 








81.0 


180.0 


90 


-81 


1.3 ± 0.7 








63.0 


180.0 


90 


-63 


4.3 ± 2.6 








9.0 


180.0 


90 


-9 


3.5 ± 1.7 
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TABLE XVII: The measured (ep — > eprf) cross section d 5 a/dk' elab d{l' elab d{l lc . m . (± statistical error ± systematic error) in 
pb/(GeV sr 2 ), for data set I-b. The (ep — > ep^) kinematics are entirely determined by the five variables: Q 2 (= 1.0 GeV 2 
fixed), beam energy E beam (= 4.032 GeV fixed), and (W, cos 6 C . m . , ip) which are given in the Table. Numbers in parenthesis 
represent a rough estimate of the systematic error, when missing. 



W (GeV) 


<p = 15° 


tp = 45° 


<p = 75° 


Lp = 105° 


Lp = 135° 


Lp = 165° 








cos c . m , = 


-0.975 






1.05 












140 ± 73 ± 62 


1.07 








69 ± 41 ± 18 


63 ± 24 ± 25 


120 ± 30 ± 24 


1.09 






50 ± 32 ± 9 


91 ± 24 ± 10 


92 ± 20 ± 12 


86 ± 19 ± 24 


1.11 




94 ± 40 ± 2 


53 ± 17 ± 9 


114 ± 20 ± 11 


110 ± 21 ± 9 


61 ± 16 ± 18 


1.13 




36 ± 14 ± 7 


54 ± 14 ± 15 


75 ± 16 ± 8 


72 ± 17 ± 17 


53 ± 16 ± 3 


1.15 


43 ± 14 ± 15 


95 ± 18 ± 15 


96 ± 17 ± 9 


97 ± 18 ± 14 


131 ± 23 ± 43 


104 ± 19 ± 27 


1.17 


55 ± 13 ± 9 


112 ± 18 ± 8 


102 ± 17 ± 12 


124 ± 18 ± 9 


179 ± 22 ± 24 


168 ± 22 ± 10 


1.19 


93 ± 16 ± 7 


116 ± 17 ± 13 


136 ± 18 ± 13 


154 ± 18 ± 23 


145 ± 19 ± 14 


178 ± 24 ± 15 


1.21 


118 ± 18 ± 17 


116 ± 17 ± 20 


167 ± 18 ± 9 


119 ± 15 ± 10 


152 ± 21 ± 6 


144 ± 24 ± 20 


1.23 


111 ± 16 ± 13 


109 ± 15 ± 18 


102 ± 13 ± 11 


141 ± 16 ± 8 


107 ± 17 ± 18 


83 ± 15 ± 14 


1.25 


51 ± 11 ± 12 


78 ± 12 ± 11 


94 ± 12 ± 11 


74 ± 13 ± 10 


81 ± 11 ± 6 


96 ± 12 ± 10 


1.27 


41 ± 9 ± 7 


51 ± 9 ± 6 


48 ± 9 ± 6 


64 ± 9 ± 4 


61 ± 8 ± 4 


47 ± 7 ± 5 








COS# c . m . = 


-0.875 






1.05 










296 ± 90 ± 18 


103 ± 37 ± 11 


1.07 








169 ± 53 ± 12 


163 ± 36 ± 6 


139 ± 30 ± 14 


1.09 






374 ± 163 ± 9 


98 ± 25 ± 40 


113 ± 25 ± 14 


122 ± 29 ± 11 


1.11 






107 ± 32 ± 8 


129 ± 24 ± 10 


119 ± 28 ± 13 


58 ± 23 ± 45 


1.13 






91 ± 21 ± 6 


123 ± 24 ± 4 


68 ± 21 ± 16 


99 ± 31 ± 11 


1.15 




81 ± 29 ± 12 


143 ± 23 ± 9 


97 ± 20 ± 18 


99 ± 22 ± 13 


117 ± 25 ± 10 


1.17 




38 ± 15 ± 15 


99 ± 18 ± 5 


106 ± 18 ± 11 


118 ± 22 ± 14 


106 ± 30 ± 14 


1.19 


172 ± 44 ± 26 


112 ± 19 ± 8 


127 ± 18 ± 11 


125 ± 18 ± 12 


145 ± 29 ± 18 


272 ± 122 ± 90 


1.21 


103 ± 22 ± 11 


123 ± 19 ± 7 


134 ± 17 ± 5 


183 ± 21 ± 11 


178 ± 51 ± 18 


217 ± 70 ± 46 


1.23 


129 ± 21 ± 8 


117 ± 18 ± 6 


135 ± 16 ± 9 


133 ± 20 ± 12 


130 ± 24 ± 42 


68 ± 19 ± 28 


1.25 


99 ± 17 ± 6 


95 ± 17 ± 6 


99 ± 13 ± 7 


68 ± 14 ± 8 


88 ± 14 ± 14 


81 ± 21 ± 8 


1.27 


54 ± 13 ± 7 


60 ± 13 ± 14 


87 ± 13 ± 5 


71 ± 11 ± 6 


57 ± 12 ± 7 










COS# c . m . = 


-0.650 






1.03 












287 ± 103 ± 11 


1.05 








201 ± 177 ± 16 


211 ± 58 ± 14 


188 ± 46 ± 9 


1.07 








191 ± 52 ± 10 


112 ± 28 ± 11 


131 ± 31 ± 20 


1.09 








150 ± 31 ± 9 


110 ± 26 ± 10 


72 ± 32 ± 101 


1.11 






114 ± 36 ± 5 


108 ± 22 ± 7 


179 ± 46 ± 41 


153 ± 77 ± (120) 


1.13 






115 ± 25 ± 5 


139 ± 27 ± 7 


62 ± 24 ± 515 


78 ± 32 ± (120) 


1.15 






86 ± 18 ± 5 


116 ± 23 ± 11 


128 ± 27 ± 169 


91 ± 36 ± (120) 


1.17 




145 ± 48 ± 13 


149 ± 21 ± 10 


137 ± 22 ± 15 


132 ± 33 ± 400 




1.19 




129 ± 30 ± 24 


169 ± 21 ± 5 


152 ± 22 ± 27 


178 ± 69 ± (200) 




1.21 




175 ± 28 ± 22 


152 ± 19 ± 14 


125 ± 21 ± 17 






1.23 


66 ± 33 ± (30) 104 ± 19 ± 16 


131 ± 17 ± 10 


87 ± 21 ± 10 


88 ± 38 ± (200) 




1.25 


79 ± 23 ± (20) 108 ± 19 ± 12 


76 ± 12 ± 12 


50 ± 15 ± 59 


125 ± 30 ± (200) 




1.27 


80 ± 19 ± (20) 107 ± 20 ± 11 


102 ± 15 ± 15 


85 ± 16 ± 73 


89 ± 42 ± (200) 





